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In this, the first of two companion papers, we present a method for finding the gravitational 
self-force in a modified radiation gauge for a particle moving on a geodesic in a Schwarzschild or 
Kerr spacetime. An extension of an earlier result by Wald is used to show the spin-weight ±2 
perturbed Weyl scalar (ipo or ipi) determines the metric perturbation outside the particle up to 
a gauge transformation and an infinitesimal change in mass and angular momentum. A Hertz 
potential is used to construct the part of the retarded metric perturbation that involves no change 
in mass or angular momentum from ?/>o in a radiation gauge. The metric perturbation is completed 
by adding changes in the mass and angular momentum of the background spacetime outside the 
radial coordinate ro of the particle in any convenient gauge. The resulting metric perturbation is 
singular only on the trajectory of the particle. A mode-sum method is then used to renormalize 
the self-force. Gralla shows that the renormalized self-force can be used to find the correction to a 
geodesic orbit in a gauge for which the leading, 0(p~ 1 ), term in the metric perturbation has spatial 
components even under a parity transformation orthogonal to the particle trajectory, and we verify 
that the metric perturbation in a radiation gauge satisfies that condition. 

We show that the singular behavior of the metric perturbation and the expression for the bare 
self- force have the same power-law behavior in L — I + 1/2 as in a Lorenz gauge (with different 
coefficients). We explicitly compute the singular Weyl scalar and its mode-sum decomposition to 
subleading order in L for a particle in circular orbit in a Schwarzschild geometry and obtain the 
renormalized field. Because the singular field can be defined as this mode sum, the coefficients of 
each angular harmonic in the sum must agree with the large L limit of the corresponding coefficients 
of the retarded field. One may therefore compute the singular field by numerically matching the 
retarded field to a power series in L. To check the accuracy of the numerical method, we analytically 
compute leading and subleading terms in the singular expansion of ipo and compare the numerical 
and analytic values of the renormalization constants, finding agreement to high precision. Details 
of the numerical computation of the perturbed metric, the self-force, and the quantity h a 0U a u 13 
(gauge invariant under helically symmetric gauge transformations) are presented for this test case 
in the companion paper. 

PACS numbers: 04.30.Db, 04.25. Nx, 04.70.Bw 



I. INTRODUCTION 



Among the most important sources for LISA are extreme mass ratio inspirals (EMRIs) of solar mass compact 
objects into supermassive black holes. LISA couldpotentially measure hundreds of EMRI events [llwhose wide range 
of astrophysical and fundamental implications 0, Q include determination of the Hubble constant [4, ,5| ; of luminosity 
distance, mass and spin of galactic black holes; and measurements of the deviation from a Kerr geometry of the central 
object [f|. Such measurements depend on accurate parameter estimation, for which it is essential to have accurate 
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gravitational waveforms available; these require an accurate calculation of the gravitational self-force experienced by 
the particle. 

The gravitational self-force contains both dissipative and conservative parts. The dissipative part is simply the 
familiar contribution from the half-retarded minus half-advanced Green's function, smooth at the position of the 
particle. The more difficult portion of calculating the gravitational self-force is the determination of the conservative 
part of the field. Estimates of its effect on the phase of the waveform are given, for example in 0-EI)- The focus 
of our work is on this latter problem. A number of authors have worked on this problem in recent years. Poisson's 
Living Review [l2j is a comprehensive self-contained introduction to the subject; Barack 13| gives an extensive recent 
review; and Detweiler L| provides a more concise summary. 

A complicating feature of the conservative part of the self-force is a combination of its inherent non-locality in 
curved space due to scatter off curvature [l5l - [l7l |) and the singular behavior of its expression near the particle. The 
calculation can be made tractable through the observation that the field consists of both a regular part that is entirely 
responsible for the self-force and a singular part that contributes nothing to the self-force. 

The other complicating feature of self-force computations is the choice of gauge. The Lorenz gauge is particularly 
useful for sorting out formal issues but does not lead to separable, decoupled equations in the Kerr spacetime. For 
this reason we choose to work in a radiation gauge that exploits the separability of the Tcukolsky equation [19, [l9j] 
for the perturbed Weyl scalar ipQ in the Kerr background. Using the Her tz p otential formalism developed for gravity 
by Kegeles and Cohen [HIHI, Chrzanowksi (H, Stewart [H and Wald [Hf (which we will call the CCK formalism), 
it is possible to construct a metric perturbation from source-free solutions to the Teukolsky equation. Explicit forms 
of the reconstruction a re g iven for perturbations of Schwarzschild and Kerr spacetimes, respectively, by Lousto and 
Whiting [2II and Ori [26j |. With ipo written as a sum of angular harmonics (in Kerr, a sum of oblate spheroidal 
harmonics), the reconstruction yields the part of the perturbed metric that has no change in the mass or angular 
momentum of the spacetime. There is a radiation-gauge form of a change in mass (an £ = perturbation in the case 
of a Schwarzschild b ackg round) that arises from a Hertz potential in the CCK formalism, but it is singular on a ray 
through the particle [27l . [HI ■ One can find a nonsingular form for the metric of a mass perturbation in a radiation 
gauge [28| . but we know of no advantage to using it. To obtain a gauge in which the perturbed metric is singular 
only on the particle's trajectory, we add the change in mass and angular in a gauge for which that part of the metric 
perturbation is continuous. 

An expression for the self-force is then found as a mode-sum in terms of the retarded metric. To renormalize 
it, one must subtract off a singular part that does not contribute to the self-force. We consider two alternatives, 
involving either an analytic or numerical determination of the singular part of the expression for the self-force f a as 
a power series in the integer i that labels the angular eigenvalues (more precisely, a series involving L := I + 1/2). 
Rcnormalization relies on subtracting from the bare self-force - from the expression for the self-force in terms of the 
retarded field - a singular part that does not contribute to the self-force. This can be checked by showing that the 
singular vector field f sa that one subtracts has vanishing angle-average over a small sphere about the particle. In 
the case of a circular orbit about a Schwarzschild black hole, the conservative part of the self force is axisymmetric 
about a radial ray through the particle. We numerically compute the axisymmetric part of / SQ and show that (as in 
a Lorenz gauge) it coincides with the axisymmetric part of — mV Q i, with p the geodesic distance to the trajectory. 
Because the angle-average of V"i over a sphere of radius p about the particle vanishes as p — > 0, the renormalized 
self-force gives the first-order correction to geodesic motion in our modified radiation gauge. 

The paper is organized as follows: In Sec. II, wc introduce the self-force equations and a criterion for their use in 
a generic gauge; we briefly review features of mode-sum renormalization in a Lorenz gauge; and we review relevant 
parts of the Newman- Penrose [19, [3(3] formalism and spin- weighted harmonics. In Sec. Ill, we begin with a list of the 
steps involved in computing the self-force in a modified radiation gauge. In the subsections that follow, we obtain a 
simple analytic expression for the singular part of each of the Weyl scalars ipo and we relate the small-distance 
behavior of the Weyl scalars to their large I behavior, and observe that the singular behavior in I of the expression 
for the self-force has the same behavior (involves the same powers of £) in a radiation gauge as in a Lorenz gauge; 
and we show that the perturbed metric obtained from ip is unique up to gauge transformations and the addition of 
metric perturbations corresponding to infinitesimal changes in mass and spin. We conclude the section by studying 
the parity of the radiation-gauge part of the perturbed metric; in particular, we show that spatial part of the metric 
(in the frame of the particle) is even under parity to leading order in the geodesic distance p to the trajectory. In Sec. 
IV, wc specialize to a particle in a circular orbit around a Schwarzschild black hole, finding tpQ Ct and ipQ, the retarded 
and singular forms of ipo and comparing the analytic and numeric methods of renormalizing ipQ. The substantial 
analytical work involved in the mode-sum expression for the leading and subleading terms of ipQ is detailed in an 
appendix. Finally, in Sec. V, we briefly discuss the results. 
Conventions 

Greek letters early in the alphabet a, (3, . . . will be abstract spacetime indices: letters //, ;/.... will be concrete 
spacetime indices, labeling components in Schwarzschild or Boyer-Lindquist coordinates. Bold-face Greek indices 
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fi, v will label components along the null Newman-Penrose (NP) tetrad defined in Eq. (jT8|) below. We adopt the 

H signature of Newman and Penrose and use the standard names l a ,n a , m a , rh a for the null NP tetrad and NP 

notation for the spin coefficients. 



II. REVIEW OF SELF-FORCE AND OF BLACK-HOLE PERTURBATIONS IN AN NP FRAMEWORK. 



A. Gravitational self-force 



We work in linear perturbation theory, for which the metric perturbation is a solution with point-particle source to 
the Einstein field equation linearized about a Kerr or Schwarzschild background. With the particle's mass, trajectory 
and velocity denoted by m, z(r) and u a , respectively, the source is given by 

T aP (x)=mu a u p f 6 4 {x,z(T))dT. (1) 



Here u a = u a (x) when x lies on the particle's trajectory, and the 5- function is normalized by / S 4 (x, z) V\g\d 4 x = 1. 

We denote by h T *g the retarded solution to the perturbed Einstein equation with this source. As noted by Quinn 
and Wald and by Detweiler and Whiting [3l| , in the MiSaTaQuWa prescription for finding the self- force [HI, , a 
particle follows a geodesic of the metric g a p + h™p, where h^a is given by 

K$ = h$-Kp> ( 2 ) 

with h B a p a locally defined singular field, chosen to cancel the singular part of h™p and to give no contribution to the 
self-force. The perturbed geodesic equation has the form 

a a := u ■ Vu a = -{g aS - u a u s ) \V^f - \v S h^ v?u>. (3) 

Here a a is the acceleration with respect to the background metric, and the self- force is, by definition, f a = ma a . 

We will denote by a rctQ the expression on the the right side of Eq. <j3j> , with h 1 ^ replaced by h™p. Work by 
Gralla [32j| . following a careful justification of the self- force equations by Gralla and Wald (33|, gives the following 
characterization of a", based on the vanishing of the angle-averaged singular part of the expression for the self- force: 
Let p be geodesic distance to the particle trajectory. Let h™p be the retarded metric perturbation in a gauge for which 
its spatial part near the trajectory is 0(p _1 ) and has even parity to that order. Then the self- force is given in local 
inertial coordinates about P by 

a rcnA1 = lim f a rcttl dn (4) 

where S p is a sphere of constant p in a geodesic surface orthogonal to the worldline. That is, the first-order perturbative 
correction to the geodesic equation is 

irv pu = a , (oj 

with a rcnQ given by Eq. (0J. 

One can thus identify the singular part of the self- force with any vector field mo s a near the particle trajectory for 
which a rct a — a s a is continuous and 



/ a s vdn = 0. (6) 
Js„ 



Then 



lim 



'(P) = ^lim p {a ret "(F') - a a [h s ](P')}. (7) 



A free particle in flat space has no self-force, and the form of its linearized field can be used to obtain a s a in a curved 
spacetime. Its linearized gravitational field in a Lorcnz gauge is the Schwarzschild solution in isotropic coordinates, 
linearized about flat space 

, 2m r 2m. . 

h^u = — Op,, = — (Tffu, - ZuLpUv). (8) 
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In a Lorenz gauge for a particle in geodesic motion in a curved background, the singular part h s aj3 of the metric 
perturbation takes this form in local inertial coordinates T, X, Y, Z centered at at any point P along the trajectory, 
with &t = u: 

h s £ 01 = -— <W = y (9^ - 2<W| P . (9) 

The corresponding singular part of the self-force is given by 

As in flat space, the angle-average of / M over a sphere of constant p vanishes in the small-p limit. In particular, 
although the leading correction to the flat-space coordinate expression x t /p 3 can be a term of order p°, the term has 
the form CjkX l x^ x k / p 3 ; because the term is odd in x % , we have 

lim [ V- dQ = 0. (11) 
p^°Js p P 

Because the self-force involves only the values of V 7 /i^™ on the particle's trajectory, two tensors h s a p and h s a p give 
the same self- force if V 7 (/i^g — h a a g) vanishes on the particle's trajectory. In particular, Detweiler and Whiting [3l[ 
show that there is a choice h^p of the singular field that is a locally defined solution to the perturbed field equations 
with the same point-particle source. One can choose local inertial coordinates (THZ, coordinates, for example) for 
which the Detweiler- Whiting singular solution differs from the form §§§ by terms of order p 2 [34|]. Following the 
notation in Detweiler- Whiting, we denote their form of the singular field by an upper-case 5". 



B. Mode-sum renormalization 

A review of mode-sum renormalization, is given in Ref. [l3| . We recall some of its main features in a Lorenz gauge; 
many of these continue to hold in our (modified) radiation gauge. 

In mode-sum renormalization, the retarded field h r °p and the corresponding expression a rctQ are written as sums 
over angular harmonics, labeled by £, m. In a Schwarzschild background, these are unambiguously associated with 
the spherical symmetry of the backgound. In Kerr, one can use the spherical coordinates of a Boyer-Lindquist chart 
to define the decomposition. When the retarded field is written as a superposition of angular harmonics, its short- 
distance singular behavior ([9]) is replaced by a large L-divergence of the mode sum at the position of the particle. 
(Appendix [B] relates the large L behavior of an function on a sphere to its singular behavior for small p.) 

For a particle at coordinate radius 7"o, the angular harmonics have finite limits as r — > ro from r < tq or r > ro- 
Denoting by h™^ the sum over m of all harmonics belonging to a given t, one has 

fCf(P) = V + B^/L + 0(L- 3 ). (12) 
Similarly, with aj et " the contribution to a rcta from h™* e , the large-L behavior of a r e eta is given by 

a Y t "(p)=A ± ^L + B' 1 + 0(L- 2 ). (13) 

Explicit expressions for the renormalization coefficients and have been found for generic orbits in a Kerr 
background by Barack who shows that the first two terms in this expansion reproduce the singular part of the 

acceleration, — V Q -, up to terms that vanish at the particle: 
P 

a s e f *(P) = A ± f l L + B> 1 . (14) 

The fact that the term of order L _1 vanishes for each component a SM is related to the behavior of a short-distance 
expansion in which terms of order p k with k even occur with an odd number of factors of x l . The retarded acceleration, 
expressed as a mode-sum that diverges on the particle trajectory, is regularized by a cutoff £ max in t, 

^max 

a rc gM = J2(af ta -a s /), (15) 

1=0 
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and the renormalized acceleration is given by 



lim Y(a? ta -a s /) 



L ax ->co 



(16) 



C. Black-hole perturbations in an NP framework 



The present method obtains the metric perturbation from components of the Weyl tensor along basis vectors of an 
NP tetrad 



e?:=i a , e£:=n Q , eg := m Q , e£ := m c 



(17) 



whose two real null vectors l a and n a are along the principle null directions of the Kerr geometry. In particular, the 
Kinncrsley tetrad has in Boyer-Lindquist coordinates the components 



m = (- 



i,o, x ), K) 



i 



(r 2 +a 2 ,-A,0,a), (m") 



1 



(ia sinf, 0, 1, ij sin#), 
(18) 



A A" v " ' 2(r 2 + a 2 cos 2 0) v v~> ^(r + iacos 

where A = r 2 — 2Mr + a 2 . We denote the associated derivative operators by 

D = l»d^ A = n"5 MJ « = m"0 M , 5 = m"0 M , (19) 
where boldface distinguishes these operators from subsequently defined scalars. The nonzero spin coefficients are 



1 



Q = 



r — iacos( 



P = -Q 



cot I 



2^2' 



7r = — j=ag smfl, r = pfflggsmt/, 

V 2 v 2 



(J-=-xQ 2 Q^, 1 = V+^QQ(r - M), a = ir-p, 



(20) 



where we distinguish g from p introduced before Eq. (U). The spin- weight s = ±2 components, ipo and "04 1 of the 
perturbed Weyl tensor are given by 



-04 = —C a j3ysn- a fh^n 1 ih s . 
Each satisfies the decoupled Tcukolsky equation appropriate to its spin weight: 



(21) 
(22) 



%ip s := 



(r 2 +a 2 ) 2 2 . 2 

— cr sur 1 

A 

1 d 



if 
or- 



- 2s 



M(r 2 - a 2 ) 



— r — ia cos ( 



AMar d 2 





- 2s 









sin 8 88 
An(r 2 + a 2 cos 2 8)T t 



i cos 8 


d 


'a 2 


1 


sin 2 8 _ 


d<j> + 


A 


sin 2 6» 



A~ s — ( A s+1 — 
dt A dtd(f> dr \ dr 

+ (s 2 cot 2 8 - s)^ s 



where 



1ps=2 = "00, 

T s=2 = 2(8 + 7f- a- 30 - At)[{D - 2e- 2g)T 13 - (S + 7f- 2a - 2/?)Tn] 
+2(D - 3e + e - Ag - g) [( 6 + 2n - 2f3)T 13 -(D-2e + 2e- g)T 33 ], 



(23) 



(24a) 



T s= -2 



4- 04, 



= 2g- 4 (A + 3 7 - 7 + Afi + p)[(8 - 2f + 2a)T 2i - (A + 2 7 - 27 + JX)T i4 \ 
+2g- i (S - f + p + 3a + 4tt)[(A + 2 7 + 2/2)T 24 -(6-f + 2p + 2a)T 22 



(24b) 
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The scalars ipQ et and g 4 ip r 4 ot can be decomposed into time and angular harmonics, 

/.ret 7-> o i(md>—ujt) /or„\ 

Q-^Tim* = -zRemu -2S imu e^- ut \ (25b) 
where the s Se mhJ are oblate spheroidal harmonics and where s Rtmu> satisfies the radial equation, 

s+1 dR\ t f K 2 — 2is(r — M)K 



A- s — \^A s+1 —j + ^ ^ '— + Aisur -\)R= -ATrT stmuJl (26) 

with K — (r 2 + a 2 )uj — ma. The source, a distribution involving 5(r — r ) and its first two derivatives, is obtained from 
the source term in Eq. (f2"3")) using the completeness and orthogonality of the spin-weighted spheroidal harmonics. 

Each angular eigenvalue A is a continuous function of a. For a Schwarzschild background, A has the value (£ — 
s)(£ + s + 1), and for Kerr it is labeled by its value of £ at a = 0. 

For perturbations of Schwarzschild, tensor components with s± indices along m a and S2 indices along fh a have 
angular behavior given by spin- weighted spherical harmonics s Ye m (9, (f) with spin- weight s = s\ — s%, where 

= f [(^- s )!/(^ + s )!] 1/2 g s r, m , o< s <£, 
s em I (-iy[(£ + sy./(£- s y.] 1/2 8- s Y tm , -£< s <o, 



with 



3?/ = — (dg + i esc Od^ — s cot 9) r?, 

h = -(dg-i esc 98^ + s cot 9) r). (28) 



D. Reconstruction of the metric perturbation from t/'o or ^4 

The CCK procedure for obtaining metric perturbations from perturbed Wcyl scalars was developed by Chrzanowski 
and by Cohen and Kegeles (20l - [22l | (see also Stewart [23l ), and a simp ler derivation is given in Wald . Discussions 
in the context of the self- force problem are found in |25l |26| and [28[ . 

The procedure gives the metric perturbation in a radiation gauge, a gauge characterized by the conditions 

l P K p = g af) h afi = 0, (29) 

or by the corresponding conditions 

n p h af) = g^Kp = 0. (30) 

Price, Shankar and Whiting [35l |36| show that a radiation gauge exists locally for vacuum perturbations of any type 
D vacuum spacetime. 

The CCK construction has two parts. Given a solution i/jq or -04 to the source-free s = ±2 Teukolsky equation, 
one first finds a Hertz potential, a function ^ that again satisfies a source-free Teukolsky equation. Then, by taking 
derivatives of the Hertz potential, one constructs a metric perturbation for which ipo a- n d ipi are the perturbed Weyl 
scalars. 

There is a striking difference between the asymptotic behavior produced by the CCK procedure and the asymptotic 
behavior of a metric perturbation in a Lorenz gauge that approximately satisfies Eq. (f29| or (|30[) . The difference is 
related to the terminology used in the literature for the two families of radiation gauges, in which the gauge satisfying 
I hap — is called the IRC or ingoing radiation gauge and the gauge satisfying n^h a p = is called the ORG or 
outgoing radiation gauge. 

Outgoing solutions in a Lorenz gauge (for example modes for which the metric perturbation has asymptotic behavior 
e -iuu j r q(j-^^ however, satisfy the IRG conditions (|2"9"|) to leading order: Because l a is along V Q u, we have 

= V h a0 = if* hap + 0(r- 2 ). (31) 

As one might expect, an asymptotically vanishing gauge vector can take one from a metric perturbation in a Lorenz 
gauge in which the IRG condition is approximately satisfied to a metric perturbation that exactly satisfies the condi- 
tion: We exhibit in Appendix [S] an explicit, asymptotically vanishing, gauge vector from a generic outgoing solution 
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h^°p in a Lorenz gauge to an asymptotically flat metric perturbation in an IRG. An analogous gauge transformation 
leads for incoming radiation to an asymptotically flat metric perturbation satisfying wrh a p = 0. 

Curiously, however, the CCK procedure yields metric perturbations in the two gauges that are asymptotically flat 
for the opposite cases: In the IRG, with l^h a p = 0, the CCK procedure yields an asymptotically flat metric for ingoing 
radiation. For outgoing radiation the CCK procedure gives a metric whose dominant components are asymptotically 
of order r. Similarly, in the ORG gauge, with n^h a f3 = 0, the CCK procedure yields an asymptotically flat metric 
for outgoing radiation. For ingoing radiation the CCK procedure gives a metric whose dominant components are 
asymptotically of order r. This then is the justification for the terms "outgoing" and "ingoing radiation gauge" 
introduced by Chrzanowski and used in the subsequent literature. 

For the gauge (|30[) . a Hertz potential ^ is related to ipo by four angular derivatives or four radial derivatives, and 
both of these alternatives are listed below. In subsequent sections, we will be concerned only with the specialization 
of these results to the Schwarzschild spacctime. In this case given a *S> that satisfies the Tcukolsky equation for i/jq, a 
metric perturbation in the radiation gauge ORG is given by 



h a p = g 4 { n a np{8 — 3a — j3 + 5ir)(5 — 4a 
~n( a rhp) [(5 — 3a + j3 + 5n + f)(A - 



4- m a mp(A + 5^t - 37 + 7) (A + fi- 47) 
47) + (A + 5/x - p - 37 - j)(S - 4a + tt)1 }* + c.c, (32) 



which we take as the starting point for the discussion that follows. Note that the sign in this equation is appropriate 

for a H signature and is opposite to that in, for example, Wald |24j |. 

In the ORG, \& is related to i/jq through four angular derivatives according to 



i> = l[£ 4 V + 12Mdt*], 

a 



(33) 



where C = 9 — ia sin 9dt- Equivalently, £ 4 = C\CqC,-\L-2 1 with C s = — [dg — s cot 9 + i esc 6d<p] — ia sin6*cV There is 
a corresponding equation involving four radial derivatives, namely 1 

1 



Q V4 



— A 2 D 4 A 2 *, 
32 



(34) 



where D is proportional to A, rcnormalizcd to make it the radial derivative along the ingoing principal null geodesies 
(the t —> —t, (f> — > —<f) version of D): 



D 



2A 



H + a 2 cos 2 



-d t + d r - —d<t 



(35) 



The corresponding equations for the (different) Hertz potential in the IRG are listed in the second line of the 
summary table below, with C := 5 + iasmOdt- 



Gauge 


Gauge conditions 


Radial equation 


Angular equation 


ORG 


n' 3 h a = O, h = 


0" 4 V>4 = — A 2 D 4 A 2 * 
y v 32 


iPo = i[£ 4 * + 12Ma t *] 
^ ■ 


IRG 


f h a p = Q, h = 


^0 = -D 4 * 


° 



TABLE I: Relations between the gauge-invariant Weyl scalars and the Hertz potentials in the two radiation gauges. 



The fact that g 4 -04 and tpo satisfy the vacuum Teukolsky equation for spin-weights T-2 when iJ/ ORG or <I< IRG satisfy 
the spin-weight ±2 Teukolsky equations follows from the relations in the table together with the commutators 



T2D 4 
71 2 A 2 £> 4 A 2 



D 4 T- 



2 

A^ 4 A 2 T 2 , T2C 4 



T-2C 4 



= £ 4 t- 
C 4 %. 



(36) 



1 D is Chandrasekhar's 3)^ and Ori's D T . Lousto and Whiting's Eq. (28) is an incorrect version of Eq. I|34|l . with A (A in their notation) 
instead of D/2. This is corrected by Whiting and Price |35| , in which A is defined as the GHP prime |37|| of D. The Ori and Lousto- 
Whiting papers also have an incorrect factor of two in each of the equations for ^ that is inherited from an error in Kegeles-Cohen 
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These are equivalent to Eqs. (40) and (56) in Chap. 8 of Chandrasekhar [38( and their adjoint relations as defined 
there. 

Explicit solutions to the equations for the Hertz potentials have been presented for a Schwarzschild background by 
Lousto and Whiting [25| and for Kerr by Ori [2(| . Ori shows that the CCK procedure gives a unique Hertz potential 
in each gauge for each angular harmonic of "J. That is, for each harmonic, there is a unique "J that satisfies both the 
angular equation and the sourcefree Teukolsky equation; there is a unique ^ that satisfies both the radial equation 
and the sourcefree Teukolsky equation; and the two solutions coincide. For ipo proportional to the harmonic 2"SWia;, 
(j/ORG j g p r0 p 0r tional to _2<%,- m ,- w - 

Note that, although Ori's metric reconstruction is done mode-by-mode, his statement of uniqueness of solutions 
does not explicitly restrict it to uniqueness of each angular harmonic. This is, however, a necessary restriction: As 
shown in Keidl et al. [28], if one requires only that ^ satisfy one of the radial or angular equations of Table iHDl 
together with the appropriate Teukolsky equation, additional freedom remains. This freedom in the Hertz potential 
corresponds to the addition to the metric perturbation of type D solutions (changes of mass and spin and addition of 
a perturbed C-metric solution) and with gauge transformations of the perturbed metric. 

With ipo and the ORG ^ decomposed in time and angular harmonics, Eq. (|33| can be inverted algebraically as 
follows, at any r outside the particle's orbit - for r m ; n < r < r max , where r m j n and r max are the perihelion and aphelion 
values of r. The harmonics of and ipo each have the form 

ipotmu = iS( m ^& l{ja ^ , (37a) 

^Imaj = 2<S , £ mw e^ m ^~ W *\ (37b) 

where R and S are solutions to the radial and angular Teukolsky equations, respectively, and R is to be determined. 
Using the identity s S[ mu} = (—l) m+s - s Si^ m ^ J , we can write the harmonic decomposition of in the form 



= (-l) m 2^-r^ -jStae""^'. (38) 
The Teukolsky-Starobinsky identity (Eqs. (9.59) and (9.61) of Ref. ^) has the form 

(39) 

where D 2 = Xfj H (X C h + 2) 2 + 8au)(m — cilo)Xch{5Xch + 6) + 48a 2 w 2 [2Ac,tf + 3(m — aoj) 2 ] and Xch, the angular 
eigenvalue used by Chandrasekhar [38[, is related to the separation constant A of Eq. (4.9) of [19[ by Xch = A + s + 2. 
Because Eq. (|33[) mixes ^ and ^ , its inversion for each angular harmonic involves a linear combination of ipotmcj and 
i-'oe -m -u> ■ We find that the algebraic inversion gives the ORG Hertz potential in the form 

*^ = 8 D 2 + 144M 2 W 2 • (40) 

We use this inversion for circular orbits in a Schwarzschild and Kerr background. 

For generic orbits, the individual harmonics ipoimuj do not satisfy the sourcefree Teukolsky equation in the region 
^min < t < fmax, where r m ; n and r max are the values of r at perihelion and aphelion; and the presence of a source 
invalidates the algebraic angular inversion. To find "J requires one to integrate one of the radial equations of Table I. 
With the Kinnersley tetrad, the IRG radial equation for '3/ has the simplest form: In Kerr coordinates u, r, 9, <f), where 
u = t — r* , with 

dr* r 2 +a 2 - f°° , 1 

and = + a dr (41) 



dr A J r A(r')' 

we have Df(u, r, 9 1 (f>) = d r f(u, r, 9, (p). The radial equation for \1/ IRG is then d^ mG ~ 2-0o> with solution 

/>OC />OG pOO pOO 

# IRG =2/ dr 4 / dr 3 / dr 2 / dnM^r,9,4>), (42) 



satisfying the vacuum Teukolsky equation for ingoing radiation when the outgoing radial null ray does not intersect 
the particle. To find at points on a t — constant surface, one can use Eq. (|4"2"j) outside the radial coordinate r of 
the particle and a corresponding integral from the horizon to r for r < ro. 
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E. Gauge transformations of the self-force 

In the form ([3]), the perturbed geodesic is parameterized so that its tangent is normalized to 1 with respect to the 
background metric. A gauge transformation of this equation was obtained by Barack and Ori (27j . and Appendix fAI 
gives an alternative, covariant derivation in terms of an infinitesimal diffco of the metric and a family of unperturbed 
geodesies. With the same normalization, changing a background geodesic by a gauge transformation generated by 
changes its tangent vector by 

S^ u a = (S% -u a up)£ € u p ; (43) 

and u a + 8^u a satisfies a geodesic equation of the metric perturbed by £^g a p. With the acceleration defined by 
6^a a := 8gurV pu a + u^V pS^u" , we have 

S^a a = -(^ - u a u p )(u ■ V) 2 ^ + fl^juV^, (44) 

and a a u a — 0. Note that the right side vanishes if £ a happens to drag a geodesic of the background spacetime to 
another geodesic of the background spacetime: This is the equation of geodesic deviation governing the connecting 
vector joining two neighboring geodesies of g a p. For general £ Q , the right side of Eq. (|44|) then measures the failure 
of an infinitesimal diffco generated by £ Q to produce a geodesic of the background metric. 

A particle in circular orbit has 4- velocity u a = u*fc Q , with k a = t a +tt<fi a a Killing vector. The perturbed spacetime 
with a particle in circular orbit is helically symmetric, symmetric with respect to k a . We show in Appendix lAI that, 
for a gauge transformation that preserves helical symmetry (for £fcf Q = 0), the gauge-transformed self- force is given 
by 

r = r + ^M^v^v^fe,). (45) 

III. METHOD FOR COMPUTING THE SELF-FORCE ON AN ORBITING MASS 

A. Overview 

In this section we outline the method for computing the self-force in a radiation gauge. Subsequent sections will 
elaborate on the details of each step of the calculation. The method is a revision of that initially suggested in Refs. [28| 
and [39j | . In broad terms, the steps involved in computing the self- force in a modified radiation gauge arc: 

A. Compute the retarded Weyl scalar, -0Q Ot (or ipT*)- 

B. Use the retarded Weyl scalar to construct a Hertz potential 4' rot as a sum of angular harmonics for r 7^ ro- 

C. Using the CCK formalism described above, reconstruct the retarded metric perturbation in a radiation gauge, 
and find the perturbation in mass and angular momentum in an arbitrary gauge. 

D. Find the expression for the self force as a mode sum involving the retarded field. 

E. Find the renormalization coefficients for the singular part f sa of the self- force and compute f rena . 

In this approach, all fields are written as a sum of time and angular harmonics. The angular harmonics (|25a[) of 
ipQ et have finite one-sided limits as r — > tq: 



1$L-= H^I™(*o,r). (46) 



r— >r. 



For a Schwarzschild background the harmonics are spin-weight 2 spherical harmonics, while for Kerr they arc oblate 
spheroidal harmonics whose form depends on the Kerr parameter a. For a Kerr as well as a Schwarzschild background, 
however, one can write the singular field in terms of spin- weighted spherical harmonics (40l - |42j . where each spin- 
weighted spheroidal harmonic is a sum of the form 

00 

s S imu e im 'f>= h'usYem, Anin = max(|s|,|m|). (47) 
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The computation of ipQ et is straightforward, involving an integration of the radial Teukolsky equation (|26|) and the 
computation of spin- weighted spherical harmonics. In computing a Hertz potential \J ,ret from "0o ct j our choice of 
radiation gauge is dictated by requiring that the CCK-constructed ^> ret vanish asymptotically: We use the ORG and 
find \l/ ret as the solution of Eq. (|4T)|) to the angular equation. The tetrad components of the metric perturbation h™% 
are then obtained from Eq. (p?2"j) and are used to compute the expression for the bare self- force in terms of h T *L 

f a [h Tet ]/m = a rcta = -{g aS - u a u s ) \VpK$ - \v S hf^ u?u\ (48) 

The remainder of the problem is the mode-sum renormalization of the self-force and the recovery of the part of the 
metric perturbation that does not arise from ipo- We discuss them in turn in the next two subsections. 



B. Mode-sum renormalization in a radiation gauge 

One can carry out a mode sum renormalization by finding the radiation-gauge version of the power series (|13|) that 
expresses the singular behavior of a ret for r > ro and r < near the position of the particle. We consider two ways 
to proceed: 

1. One can find an analytic expression for a sa as a sum in powers of L, starting from an expression we derive below 
for ipQ, the singular part of i/jq^- This analytic way follows the steps just listed in describing the path from ^Jj et to 
a rotQ , to successively obtain expressions for \I/ S , h s aS and a sa . 

2. The second method is significantly simpler: One simply numerically matches a power series in L to a^ et ' M (P) for 
successive values of I. 

We begin with a discussion of the analytic method. We find an explicit expression for -0q to sub-leading order 
in the distance to the particle's trajectory in terms of components of the tetrad vectors along and orthogonal to 
the trajectory. We then characterize the powers of L that appear in the power series for a rctsa . Finally, we turn 
to a description of renormalization by numerical matching. In Sec. IV below, we test the numerical method in a 
relatively simple case by comparing analytic and numerical renormalizations of ipo for a particle in circular orbit in a 
Schwarzschild background. 

Note that, in a mode-sum renormalization, the singular parts of the the perturbation in the metric and the self-force 
are determined by the large- L behavior of the retarded fields. They are, in particular, independent of the choice of 
gauge in which one describes perturbations of mass and angular momentum. 



1. Analytic method 

The analytic method begins by finding an analytic expression for ipQ or The decomposition of the metric 
perturbation h™g in a Lorenz gauge, 

7 ret, Lor jren,Lor . j s.Lor /,in\ 

Kp =Kp + h af s . ( 49 ) 

gives a corresponding gauge- invariant decomposition of the perturbed Weyl scalars. From the expression for the 
perturbed Weyl (or Riemann) tensor in terms of the perturbed metric, we have 

=Otfh%, % =Ofh B api (50a) 

=Ofh$, V! = Of h s a p, (50b) 

where 

Of = ^(m a m l3 Pl s + l a l f3 m~ f m s - ^m^m^l 5 - m a l^Pm s ) \7 y \7s 7 (51a) 

Of = - (fh a m^n / n s + n a n l3 rn<m s - n a fh^rh^n 5 - fh a n' 3 n' 1 fh 5 )\7 1 \7s. (51b) 

Although we will not need the Detweiler- Whiting form of the singular field, it is worth noting that using it would 
yield a smooth version of tpQ Cn (or ■04 en ) that satisfies the sourccfrec Teukolsky equation in a neighborhood of the 
particle's trajectory. That is, because h^g satisfies the perturbed field equation with the same source as h™g, the 
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corresponding field -0q is a local solution to the s = 2 Teukolsky equation with the same source as ip™*, implying 
that -0g et — -0g satisfies the corresponding homogeneous equation. 

Denote by e := \JT 2 + X 2 + Y 2 + Z 2 the distance with respect to the Euclidean metric dT 2 + dX 2 + dY 2 + dZ 2 
of the local inertial coordinates introduced before Eq. ((9]). We will argue below that the self- force in our modified 
radiation gauge, like that in a Lorenz gauge, has dominant singular behavior of order e~ 2 and can be regularized by 
subtracting leading and subleading terms in e. These arise from the leading and subleading terms in e in ip^. Instead 
of directly using Eq. (|50a[) to compute ipQ, it is easier simply to observe that Oq V and h^ lL , have to subleading order 
in e their flat space form. It follows that the value of has at subleading order in e its form for a perturbation of 
fiat space. That is, in terms of T and p, ipg is the linearized Schwarzschild field of a static particle: 

C U lS = ~y (*M + 3<^V^V 5 V - 3$V ffl T V^T - 6V [Q p V 0] T V^p V 5 ^) . (52) 
From this we obtain 

^ = _^nT m p_ lPm T^ 2 ^ (53) 

pi 

r = _^( n T m p_„p m rxa (54) 

where l T := l a V a T, l p := PV Q p, m T := m Q V Q T, m p := m a V a p. One can, of course, also obtain ([53]) and ([Ml) by 
applying the operators O^f and Cg' 3 of Eqs. (|50a[) and (|50b|) to the singular metric ©. Eqs. (|53p and (|5^|) are valid 
in general Petrov type D spacetimes. 

For each of the quantities ipo, ^P, h a p, and a", we will use the subscript I to denote the sum over m of the contributions 
from angular harmonics associated with 1, to, and we suppress the index oj. The equations that relate ipo and ^ do not 
mix spin- weighted spheroidal harmonics, and ^ B is most simply found as a sum of these harmonics. The equations, (|32j) 
and ([3j, that relate ^ to h s a/3 and /i^ to a a , however, mix spheroidal harmonics on a Kerr background. The large-L 
expansions of their components can be found in terms of spin- weighted spheroidal harmonics, spin- weighted spherical 
harmonics, or spherical harmonics. Different choices lead to different definitions of the subleading contributions, 
because of the mixing of different values of I in relating, for example, spin-weighted spherical harmonics to spin- 
weighted spheroidal harmonics. What matters is only that the same convention is used for the angular harmonics of 
the retarded and singular fields. 

Because tJjq involves two derivatives of h 8 a p and can be computed from the Lorenz singular field, its large-L behavior 

is greater by two powers of L than the large-L behavior at r = r of h s ^ m of Eq. ([T2")) : 

%f ~ E $Mm(*o) 2Y em (0o,M = A±L 2 + B±L + 0(L- V ). (55) 

m 

In Sec. IV and Appendix [Cl we find the large-L expansion of ipo for a particle in circular orbit in a Schwarzschild 
background (restricting consideration to the part of i/jo axisymmetric about the position of the particle). 

Because the Hertz potential ^ rct involves four integrals of ipo, its singular part has to subleading order in L the 
form 

*J = A±/L 2 + B±/L 3 . (56) 

This behavior also follows from the explicit form (j40|) of ^i mU j, together with the fact that, for large £, the spheroidal 
eigenvalues A approach their spherical values, X/£ 2 — > 1. The leading term in the expansion is then immediate from 
the leading term in (|55p . while subleading terms involve the expansion of A (found analytically or numerically) in 
terms of L. 

The metric perturbation involves two derivatives of \E', implying that the singular and retarded fields in a radiation 
gauge have the same leading power of L as their Lorenz counterparts, 

Finally, the self-force J rcn > R Ga j g computed from V 7 /i 1 Q ™' RG , using Eq. ([3J. As in the Lorenz gauge, the additional 
derivative gives the behavior 



gSRGfi = A y.± L + B u,± + (L- 2 ). 



(58) 
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The renormalized acceleration at the position of the particle is then given by 

a renRG M= y (ret RG M _ fl s RG M j_ (5Q) 

As noted in the introduction, for a particle in circular orbit in a Schwarzschild background, we find that the large 
L expansion of a sRGa agrees through O(L ) with 

a sRGa = (60) 
P 

differing only by a constant from its form in a Lorenz gauge. This form of a sRGa does not contribute to the self- force, 
because Eq. © is satisfied - the small- ,o angle average of a sRGa vanishes. Thus in this case, we can identify the 
singular field with its leading and subleading terms as a power series in L, 

a sRG M = A M±£ + £M±. (61) 

For a CCK radiation gauge, there is as yet no general proof that a sRGa is given by its leading and subleading terms. 
This is true in a Lorenz gauge, and we expect it to hold for a radiation gauge as well, with an argument based on 
the common property that a sa can be expressed (for r > ro or r < ro) as a power series that begins at 0(e -2 ) and 
involves positive powers of the coordinate differences x M — Xq and odd powers of p. 

Without a proof, one must check that the computed a sa satisfies Eq. J6]). 

In the Schwarzschild example below, the coordinate expression for to subleading order is a sum of more than 25 
terms, and to find its large-L expansion we computed the large-L expansion of each term. Finding the corresponding 
large-L expansion for a s RG/i involves finding a large-L expansions of all combinations of three derivatives of each of 
these terms for each component of a sRG/x . Without significant insight, this would mean finding the large-L expansion 
of about 300 terms. Given the simple form that a sa takes for the Schwarzschild circular-orbit there may be a similar 
form in the generic case and a much simpler way to find it. In its absence, the numerical method is much easier to 
use. 



2. Numerical method 

The renormalization coefficients occuring in Eq. (|58| for a s e RGai arc coefficients in the large-L expansion of a rot RGai 
evaluated at r = . Consequently, once one has found the numerical values of a ret RG/1 , it is not in principle necessary 
to carry out an additional analytic computation of a^ RGM . Instead, one can match to the sequence of values a^ etRGA1 
a power series in L of the form 

n F fj, 

a f* = A» ± L + B» ± + J2 jk> ( 62 ) 

k— 2,evcn 

finding the E£ that yield a best fit. A numerical check is that subtraction of the order L~ k term reduces the order 
of the series by one power of L. And one should check that the reduction of order holds for values of L larger than 
those used to obtain the coefficients in the matching. Finally, one checks numerically that f TC s^ converges to a value 
jren./j, as the cutoff £ max increases. The disadvantage of the numerical matching method is that, for a given desired 
accuracy in a rena , one must compute a retQ to higher values of I than is required when one or more renormalization 
coefficients are known analytically. 

To identify a\ a with J 4 /i± L + B^, one must again show that a sa satisfies Eq. ©. This can done if one can find 
for each component a M expressions in terms of p and the coordinate differences — Xq whose large-L expansion is 
A^L + B^, and if this expression satisfies (|6]). If the limiting angle- average has a finite value, that finite value must 
be added back, in accordance with Eq. @, to find the self- force. 

Equivalcntly, if the first two terms in this L-expansion correspond to the terms of order e~ 2 and e" 1 in the position- 
space expansion, so that the terms involving correspond to terms of order e and higher, then these latter terms 
sum to zero. We suspect this is the case for a radiation gauge, because the behavior of the singular field as a power 
series in L implies behavior in e that is no more singular than in a Lorenz gauge and corresponds for r > tq and 
r < ro to a power series in e. This is consistent with our numerical construction of the singular part of the self-force 
for a particle in circular orbit in Schwarzschild. 
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With an analytic knowledge of the first term or terms in the expansion, the numerical method is still useful in finding 
subsequent terms, and this has been done to speed convergence in the self-force computations involving mode-sum 
renormalization in a Lorenz gauge. In particular, once one knows that can be identified with the leading and 
subleading terms A M± L + B' 1± , the computation of is given by 

final / n \ n oo 

a ren M= ^ a ret M _ A M± L _ B M±_ ^ ft + ^ E g ^ L~ k , (63) 
t=Q \ fc=2,cvcn / fc=2,even £=0 

with an error of order £max • 

To justify the numerical method, we present in Sec. IIVI below a comparison of the analytic and numerical de- 
termination of renormalization coefficients for the axisymmetric part of tpo and their contribution to the self-force. 
The companion paper presents the full numerical computation of the self-force for a particle in circular orbit in a 
Schwarzschild background. Checks of the work include a numerical computation of a quantity h a pu a u^ that is invari- 
ant under helically symmetric gauge transformations and has previously been computed in Lorenz and Regge- Wheeler 
gauges. 



C. Mass and Spin: The Remaining Metric 



The CCK reconstruction of the metric perturbation from ipo gives a perturbed metric for which there is no change 
in the mass and angular momentum. For a Schwarzschild background, this is immediate from the fact that ?pQ et 
involves only values of I with I > 2, because the construction of the perturbed metric preserves the value of I. For a 
Kerr background, both ipQ Ct and \D' ret are a sum of spheroidal harmonics with spin-weight 2, and their expression (|47[) 
in terms of spin- weighted spherical harmonics involves only harmonics with I > 2. A mass perturbation requires an 
£ = part of the asymptotic perturbed metric at 0(r~ r ). The operator in Eq. (|32|) mixes different values of t, but it 
differs from its Schwarzschild form by terms smaller by 0(r _1 ) than its leading terms. That is, both ipQ 0t and ^ ret 
are 0(r -5 ), the spherically symmetric part of the operator (|32p is 0(r 4 ), and terms that mix different values of I are 
0(r 3 ), implying h™g' has no £ = part at 0(r _1 ). 

A perturbation of angular momentum requires an I — 1 contribution at 0(r~ 2 ) in /113. To see that it also vanishes, 
we show that the correction in a Kerr background to the Schwarzschild expression for h™£ is asymptotically 0(r~ 3 ). 
The operator that relates to * in Eq. {32]) is (S-3a + ^ + 5n + f)(A + n-4'y) + (A + 5fj,-Ji-3^-^)(S-4a + n). 
Now D is a derivative along an outgoing null ray. In the Kerr coordinates of Eq. (|41[) . Df(u, r, 8, <f>) = d r f(u, r, 9, </>). 
Because each time harmonic of $ has the form $ = S(9, 4>)e~ iLJU r~ 5 {l + 0(r -1 )], D$ = S(9, (/>)e -4Wtl 0(Y~ 6 ); because 
the spin coefficients arc 0{r~ r ) or smaller, h r ^ falls off like r -2 , and the correction to its Schwarzschild behavior is 
at 0(r- 3 ). 

To complete the metric reconstruction one needs to add the contributions from a change in the mass and angular 
momentum of the spacetime outside r = r . Satisfying the perturbed field equation also requires a discontinuous 
gauge transformation associated with a change in the system's center of mass. That nothing further is needed follows 
from an minor extension of a theorem by Wald that implies that a perturbed vacuum metric is determined up to 
gauge transformations and the addition of a Petrov type D perturbation of the black- hole geometry 43j . There are 
four kinds: (i) an infinitesimal change in the black hole's mass and (ii) in its spin; (iii) the perturbative version of the 
C-metric and (iv) the perturbative version of the Kerr-NUT solution. The type D perturbations are all stationary 
and axisymmetric, and only the mass and angular momentum perturbations are smooth in the region exterior to the 
black hole: The C and Kerr-NUT perturbations are each singular on their axis of symmetry, coinciding for a Kerr 
background, with the axis of symmetry of the Kerr geometry. If we were dealing with a source-free perturbation, 
regularity at the horizon and at infinity would rule out the addition of Kerr-NUT and C-metric perturbations, and, 
after a choice of gauge, we would be left with changes in mass and angular momentum (and gauge transformations). 
To extend the argument to in our case, we note that smoothness of each time-harmonic of ifjQ 0t for r 7^ ro, 
together with the explicit form (|4T))) . implies that \& ret is smooth for r / ro. That is, smoothness of -0o ct implies that 
the coefficients of each angular harmonic in its decomposition fall off faster than any power of L. Eq. (|40[) implies 
that the coefficients of the angular harmonics of ^ fall off still faster in L. Thus each time harmonic of VP is smooth 
and has no contribution from a C or Kerr-NUT perturbation. 

The contributions to the retarded field of mass and spin were examined by Detweiler and Poisson |4|| and by 
Price [45|]. With the Hertz potential restricted to t > 2 (for a Schwarzschild or Kerr background), there is no local 
contribution to mass and angular momentum. For r < ro, this is the appropriate solution for the retarded field. 
For r > ro, one has to determine in any gauge four parameters corresponding to changes in mass and spin (two 
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correspond to a change in the spin direction and arc gauge transformations); and three parameters associated with 
gauge transformations for r > r that eliminate an asymptotic dipole. The change in mass and in the magnitude of 
angular momentum along its direction in the background Kerr spacetime can be found from the integrals 

SM = ( {2T a p - &%T)lPdS a , (64a) 
Js 

SJ = - I T a ^dS a (64b) 
Js 

over any hypersurface intersecting the trajectory. (For a Schwarzschild background, all components of S3 are de- 
termined in this way.) The remaining parameters are determined by jump conditions in the field equations across 
r = r . 

Although the metric perturbations corresponding to changes in mass and spin can be written in a radiation gauge, 
we do not see a good reason to do so. In particular, the radiation gauge form of a mass perturbation that arises from 
a Hertz potential has a singularity on the radial ray through the particle. (There is an alternative radiation-gauge 
form of a mass perturbation that is nonsingular on the axis of symmetry (39j , but it has no obvious advantage over a 
mass perturbation in another gauge.) Instead we use for h™p a radiation gauge for £ > 2, together with an arbitrary 
convenient gauge for £ < 2. 



D. Leading order parity of ipo, and h a p 

We now consider the parity of the perturbed metric in a radiation gauge. Gralla's criterion is that the projection 
of the perturbed metric to a surface orthogonal to the particle's 4-velocity is even under parity to leading order 
in p. Parity here means the locally defined diffeo that exchanges points at the same geodesic distance on opposite 
sides of each geodesic orthogonal to the particle trajectory; and the associated hypersurface spanned by geodesies 
orthogonal to a point P of the trajectory is the surface onto which h a p is projected. Because the invariance under 
parity is only to leading order, we can define parity in terms of local inertial coordinates T, X, Y, Z as the diffeo 
V : (T, X, Y, Z) i — y (T, —X, — Y, —Z). The Weyl tensor has to leading order in p its flat-space form, implying that it 
is even (invariant) under both parity and time-reversal, where time reversal, T, exchanges points with coordinates 
±T, X, Y, Z. (Its behavior under both parity and time-reversal will be needed in the argument below.) 

We first show that these symmetries are retained by the singular form of ipo, given by Eq. (|55|) . Because the tetrad 
vector fields l a and m a (and n a ) are smooth, they are constant near a point P on the trajectory to leading order in 
p (as long as the particle is not on the 8 = 0, ir axis, where m a is not defined). The corresponding scalars ! Q V a T and 
to" V a T are then even to leading order under parity. Because the cartesian components of V a p have opposite signs at 
diametrically opposite points (T, X, Y, Z) and (T, —X, —Y, —Z), the scalars l a V a p and m a V a p are odd under parity. 
Then (l T m p - lPm T ) 2 is even and hence ipo is even under V to leading order in p. 

Similarly, the scalars l a V a T and m Q V Q T are odd to leading order under time-reversal; and PV Q p and m a 'V a p are 
even, implying that (l T m p — l p m T ) 2 and ipo are even to leading order under time reversal. 

To show that ^ is even under parity to leading order in p requires additional steps. First, because 4" is obtained 
from ipo as a sum of angular harmonics, we use the fact that the leading order parts of ipo and \t in p for small 
p are associated by the transform to angular harmonics with the leading order in L part of its angular harmonics, 
for large L, as described in Appendix [Bl In particular the leading terms in p of ipo and are, respectively, 0(p~ 3 ) 
and O(p), and they correspond to the large-L terms of order 0(L 2 ) and 0{L~ 2 ) in the angular harmonics. Second, 
Eq. (|4"0|) , expressing 'J in terms of ipo, involves angular harmonics of ipo on a surface of constant t, a surface that is 
not perpendicular to the particle trajectory; because of this, the plane tangent to the surface is not invariant under 
parity. It is, however, invariant under VT, implying that the restriction of ipo to a constant t surface is invariant to 
leading order under VT (Restricted to the constant t surface, VT is a parity transformation about P.) It follows that 
^ s and h s a/3 are also invariant to leading order under VT- We give the argument in terms of the the angular equation 
(|33|) for The right side of this equation is dominated for large L by its first term, having to leading order the form 

iPo = ^ 4 *; (65) 

that is, because C 4 is, for large L, quartic in A and oj, it dominates the second term. (One cannot look at the large L 
limit with lo fixed, because uj is not independent of L: For a circular orbit, for example, lu = mO.) For the particle 
at do 7^ 0,7r, C has to leading order in p the form C = do + icsc6od$ — iasin^o^t- in Boyer-Lindquist coordinates, 
VT is given to leading order in p by (t, r, 0, <p) — > (2t$ — t, 2ro — r, 29q — 6, 2<po ~ 0), implying that each term in this 
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leading form of C is odd under VT . Then £ 4 is even to leading order: 

C 4 VT = VTC 4 . (66) 
The parts of ^ that are even and odd under VT then have sources that differ by one order in L: 

i(l + PT)[£ 4 (# odd + * cvcn ) + UMdtV} = £ 4 § cvon [l + 0{L~ X )] = 8^o von , (67a) 
^(l-VT)[£ A ($ odd + V even ) + 12Md t y] = £ 4 § odd + 0(§ cvcn x L 3 ) = 8i/jQ dd . (67b) 

With a source smaller by one order in L, the algebraic inversion (|4T))) then gives <]/ odd smaller than ij/ even by one power 
of L. 

Next note that, because p is independent of T and the tetrad vectors are smooth (and hence constant to lowest 
order in p), tpo is independent of T to lowest order in p. That is, translating i/jq by AT changes it only by a term of 
order i(jqAT, from the O(AT) change in the tetrad vectors. Now time-translating ^ from a point on the t = surface 
to a point on the relatively boosted T = surface through the same point P of the trajectory involves a translation 
by a time AT proportional to p and hence changes \& only to sublcading order in p. Thus ^ is invariant under V to 
leading order in p. 

Finally, a tensor T a p is invariant under if T a p = V*T a p, where the pullback V*T a $ has components in coordi- 
nates {x»} given by {V*T)^(Q) := d il V a d v V r T ar [V{Q)\. (We have used V^ 1 = V .) In terms of the coordinates 
(T, X, Y, Z), the requirement that the spatial projection of h a p is invariant under parity is equivalent to the condition 
that each spatial component hij is even under parity 

hij(T,X,Y,Z) = {Vh)ij(T,X,Y,Z) = (T, —X, —Y, —Z). (68) 

That this condition is satisfied follows from Eq. ([32]) for h a p in terms of ^ and the fact that the leading, 0(p _1 ), part 
of h a p comes entirely from terms quadratic in the derivative operators, in A, 8 and 5. Because ^ is independent of 
T to leading order in p, each derivative operator along a tetrad vector involves only the spatial (X 1 ) components of 
the vector, implying that the quadratic derivatives A 2 ^, . . . , S 2, $> all have even parity to leading order. Finally, the 
lowest-order constancy of the tetrad vectors implies that the products of their components riiUj, . . . , m^fhj) in local 

ret RG 

inertial coordinates are even to leading order in p under parity. We conclude that the projection of ' orthogonal 
to the 4- velocity is even under parity at leading order in p: to 0(p~ 1 ). 2 



IV. PARTICLE IN CIRCULAR ORBIT IN A SCHWARZSCHILD GEOMETRY 



As a simplest explicit example of the method, we consider a particle of mass m in circular orbit at radial coordinate 
ro about a Schwarzschild black hole of mass M . In this section we compare the numerical and analytic renormalization 
methods by looking at the mode-sum renormalization of the axisymmetric part of ipo. We first compute the retarded 
field; we then find an analytic expression for the singular field to subleading order; finally, we obtain the renormalization 
coefficients of the singular field numerically, finding agreement to high accuracy with their analytic values. The 
numerical computation of the self-force is described in the companion paper. 

We work in Schwarzschild coordinates and adopt the notation 

ds 2 = fdt 2 - r X dr 2 - r 2 (d9 2 + sin 2 6d(j) 2 ), (69) 

where 

/(r):=l-— • (70) 
r 

The Kinncrsley tetrad vectors have components 

(H = (V/(r), 1,0,0), K) = i(l,-/(r),0,0), (m») = -L(°> °> */ sin ^- ( 71 ) 

* V 2r 



2 In fact, the argument shows that each component h^ v , regarded as a scalar, has even parity at leading order in p. 
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With this choice of tetrad the nonzero spin coefficients are 

Q=-~, P = 
r 

with corresponding Christoffel symbols 





cot V 




71 r 

M 


1 


—a - 




7 = 




Li 


1 12 


- -r 2 

— 1 22 


= 2 7 






1 43 


— 1 34 — 


r 4 - 

1 24 — 






r 3 

1 13 


_ p2 _ 
— 1 43 — 


r 2 - 

1 34 — 


r i4 = ~Q 




r 3 

1 33 


_ p4 _ 
— 1 44 — 


1 43 


= -rL = 2/3. 





2M 
r 



The only nonzero component of the background Weyl tensor is 

M 

* 2 = t • 



The particle's 4- velocity is 



(72) 



(73a) 
(73b) 
(73c) 
(73d) 



(74) 
(75) 



with t a and </>" timelike and rotational Killing vectors and with u* = yl — SAI/tq. Its energy and angular momentum 
per unit mass, _B := —u a t a and J := u a (f> a , are given by 



£ = 



r - 2M 



J = 



V^g - 3Mr r Q - MI 

From Eq. (|32|l . the nonzero components of the metric perturbation are 



(76) 



/in 



/133 = -r 4 



-(8 2 * + S 2 *), 

"1 (a? - w + / 2 « - + + r2 " 2M2 



2r 2 



2r 2 



To compute the self-acceleration ([3]) in terms of these tetrad components, we use the relations 



t a = -fl a + n a , 
2 



- — IT 

( m <*-fh a ), V a r = --fl a + n a 1 



to write 



(77) 
(78) 

(79) 
(80) 



Then, expanding the covariant derivatives and using Eqs. (|72|l and (|73p . we find 



(81) 



« r = («*) a (/o 2 

/o 



16- 70 8 8 v y 2r 2 

1 M _ 1 M 1M\, 

- — £> - 7 — rA + -— h 33 + 

8 r 4 r /o 2 r 2 / 



7I^ A -47I^ (s - §) + 271 fi 



.(82) 



A. The retarded field 



The retarded fields ifjf and </'y ct are simplest to compute in coordinates for which the unperturbed orbit lies in the 
6 = ir/2 plane. The particle's trajectory is then given by <f> = fii, where ft = yM/rg. From Eq. ([1]), its stress-energy 
tensor is 

m 



T Q/3 = -l- s u a u p 6{r-ro)S(coB{0))S(4>-SU) 



V ^u a u p 8(r - ro) s Y im {6, 4>) s Y lm (J, fit) . 



(83) 
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Because the source and the background metric are both helically symmetric, Lie derived by t a + fl(f> a , the retarded 
fields - metric and Weyl tensor - will also be helically symmetric. Because the tetrad vectors (|7Tj) are also helically 
symmetric, the symmetry is shared by the scalars ifjQ Bt and V'FS which therefore only involve <j> and t in the combination 
4> — fit. In the harmonic decomposition, <fi and t then occur only in the combination e lm (0- f2t ) ] an( j the frequency 
associated with each value of m is uj = mfl. 

The scalars ^l et and ipo et satisfy the Bardeen-Press equation [Z(| , the a = form of the Teukolsky equation, namely 



T s tp 



A * 



Mr 2 \ A _ s d 



A 



dr 



ip = 4irr T s 



We will work with tpo et ' whose source, from Eq. (|24a|) . has the form 

T s=2 = -2(6 - 2p)8T xx + 4(£> - Ag)(S - 2/3)T 13 - 2{D - 5g)(D - g)T 33 
=: T (0) + + T (2) 



(84) 



(85) 



where the superscripts indicate the maximum number of radial derivatives in each term. From Eq. (|83[) , these terms 
have the explicit forms 



T (0) 
T (i) 

T (2) 



- r )[(£ ~ W + m + 2)] 1/2 2Y em (9, 4>)Y im (J, 



fit 



e.m 



£mfiv 



iS'(r — r ) + 



S"(r-r ) 



mfl 4i\ 
fo r J 



5(r - r ) 



[(£ -!)(£ + 2)}^ 2 2 Y lm (9, ^tYtm (|, fit 



(86) 
(87) 



6 2imfl\ , 



ro 



fo J 



S'(r - r ) 



m 2 fl 2 6imfl 2imMfl 4 . 
—ff~ + -^) s ( r - r o) 



{if A 



)2Yim ( 



ro/o r 2 f 2 r 2 , 
Each mode of ipo or ip^, 

i> = e- iut R(r)y em (e,<t>), 
has radial eigenfunction R or R4 satisfying the radial equation corresponding to its spin- weight: 

~u 2 r 4 4iujr 2 (r - 3M) 



AR' ' + 6(r - M)R' 
AR'l -2{r-M)R' A 



A A 

W 2 r 4 4iwr 2 (r - 3M) 



(£-2)(i + 3) 
(£ -!)(£+ 2) 



Ro = 0, 
i? 4 = 0, 



(89) 

(90) 
(91) 



where the prime denotes a derivative with respect to the radial coordinate r. Solutions to these equations are related 
by 



Ro 



R-i 



To compute ipo et ' ^ is helpful to define a Green's function G(r,r') as the solution to 



AG" + 6(r — M)G' + 



wV 4iwr 2 (r-3M) \ - S(r - r^A 1 / 2 



A 



(*-2)(* + 3) G 



namely 



G(r, O = " E ^™ LA ^^<)^oo(r>), 



(92) 



(93) 



(94) 



where i?H and i?oo are solutions to the radial equation for -0o that are regular at the horizon and at infinity, respectively, 
and the quantity 



1 



Urn 



A3(i? H i?^ - i?ooi? H ) 



(95) 
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is a constant, independent of r. The full spatial Green's function G(x, x') = G(r, 9, <j>; r', 9' , <f>') is then given by 

G(x,x') = -Y. A ^^ r jy RH(r < )R oo (r > ) 2 Y em (0,<l>)2Yi m (0',<l>'). (96) 

Cm 

The Weyl scalar ip is defined and smooth everywhere except on the trajectory of the particle. It is given in terms 
of the source and the Green's function G(x,x') by 

0o = 4tt / T(t,x')G(x,x')r' 2 dV 



An 



(Ti 0) +T^)G{x,x'y 2 dV' 



=: 4 0) +4 1] +ri 2 \ (97) 

where the superscripts on tpo correspond to the three terms in Teukolsky source function defined in Eq. (|85|) . The 
three terms in Eq. (p?Tf have outside the particle trajectory 3 the form 

<> = 4W^^A, m [(^-l)^+l)(£ + 2)]V 2 i? H ( r< )^ (98) 

tm 

4 X) = 8nimnu t A Y,A im [(£-l)(£ + 2)} 1 /\Y em (e,^) 1 Y em {^,n?j x (99) 

Cm 

{[imilr 2 + 2r }R H (r < )R oo {r > ) + A [i^(r )iioo(r)fl(r - r ) 
+i?. H (r)i?^(r )%o-r)]}, 

4 2) = -47rmOV^^ ro2 ^ ro (e,0) 2 y« ro (| ) Oi) x (100) 

[30r^ - 80Mr 3 + 48Af 2 rg _ m 2^2 r 6 _ 2A 2 _ 24A ro(r - M) + 6imQr^{r - M)]i? H (r < )i? 00 (r>) 
+2(6r 5 - 20Mr 4 + 16M 2 r 3 - 3r A 2 + imQA rt)[RH(ro)Roo(r)9(r - r ) + R^(ro)Ra(r)e{ro - r)\ 
+r 2 A 2 [R^(r )R oo (r)9(r - r„) + R^(r )R B (r)9(r - r)]\ . 



B. The singular field 

Because the conservative part of the self-force is radial, it is axisymmetric about a radial ray through the particle. 
We will compare the analytic to the numerical determination of the singular part of a Weyl scalar by looking at the 
axisymmetric part of tpQ and (in Appendix [D]) its contribution to the self-force. We outline the calculation of the 
leading and subleading terms in the axisymmetric part of the singular field ipQ as a sum of angular harmonics 2^0 (©) 
whose coefficients are polynomials in £, with angular coordinates and $ chosen so that the = line (at fixed t) is 
the radial line through the particle. Details of the conversion from a small distance expansion to a large L expansion 
are left to Appendix [Cj 

The analytic expression for the resulting renormalization coefficients is then compared to a numerical determination 
by matching the retarded field to a power series in L. Remarkably, although the subleading part of ipo is a lengthy 
expression - Eq. (|C1|) . we will see that its axisymmetric part, written as a sum over angular harmonics, vanishes. 
Because the angular harmonics ^Xgm are complete in Li(S 2 \ this means that, as a distribution, ?/>q has support at 
6 = 0, where 2 Y m (®) = 0. 

The expression for the retarded field is simplest for coordinates in which the orbit is in the 9 = it/2 plane. Expressing 
the singular field of Eq. (|53p as a sum of angular harmonics is simplest if angular coordinates and $ are chosen 



3 Note that the formal integral of the Green's function also gives a <5-function contribution with support on the trajectory, namely 
-47rmf2 2 u t /o~ 1 ' 5 ( ? ' - r )5(cos 0)5(0 - Qt). 
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with the particle at 9 = 0, as we have just described. To compute the difference ipQ Bn = -0g et — ipQ, one must rotate 
?/>o et to the coordinate position of or tpQ or vice-versa. Following the conventions of Detweiler et al. [13] (henceforth 
DMW) , 9 and $ are related to 8 and ^ by a rotation of the form 

sin 6 cos <j> = cos 9, 
sin 9 sin <f> = sin 9 cos <!>, 

cos 6 = sin6 sin (101) 

With the usual association of Cartesian coordinates x, y, z to r,6,<f> and of x,y,z to r, 9,$, the map is x = z, y = 
x,z = y. 

Eq. ([53)) for iffi involves the components l T = l a V a T and l p = l a V a p. We obtain these to sublcading order in 
terms of Schwarzschild coordinates: That is, with e the distance from the particle's position P with respect to the 
positive-definite metric g a p + 2u a up, the leading and sublcading terms in T and p are 0(e) and 0(e 2 ), respectively. 
The corresponding leading and subleading terms of l T and l p are then 0(1) and 0(e). 

Expansions of p and of the local inertial coordinates T, X, Y, Z about a point P in terms of Schwarzschild coordinates 
are given, for example, in Ref. |47j |. To subleading order, T has the form 

T = (E(t-t ) - Jsin9cos$) + [ (t - t )(r - r ) - — (r - r )sin6cos$ | . (102) 

Vofo r J 

It is convenient to work with p 2 instead of p; to subleading order, we have p 2 = p^ + p^ z \ where the order e 2 and e 3 
contributions are, respectively, 

(2) = ( r ~ r o) 2 + ( 2 + j2 CQS 2 gx gin 2 q _ 2EJsin 6 cos <f>(t - t ) + ^-h(t - t ) 2 , (103) 
fo r$ 



f11 M ( 2J 2 \ , , 2 . N 2JE(M-r ), w , ^ M , 

P (3) = — 1 + — (t-io) 2 (r-r ) + f 2 (*-*o)(r-r ) sin 9 cos ^^(r-r ) 3 

2^? 2 r 

+r sin 2 0(r - r ) + 2r sin 2 9 sin 2 $(r - r ) H 7— ^ sin 2 9 cos 2 $(r - r ). (104) 

Jo 

We use Eq. This expression omits <5- functions with support at the position of the particle. These do not 

contribute to the renormalization of the retarded field if the rcnormalization is done by subtracting the singular 
field from the retarded field in a neighborhood of the particle, averaging over a sphere surroun ding the particle, and 
then taking the limit as the radius of the sphere shrinks to zero (the Quinn-Wald prescription [16()- In a mode-sum 
regularization, we discard 5-functions with support at the particle in both the singular and the retarded field. 

The background Kinnersley tetrad written in terms of 9 and $ is 

j, 1,0,0^, m a = -L ^0,0,1, ^) , where again /= M - (105) 

We now expand the needed tetrad components to subleading (quadratic) order in Schwarzschild coordinates, about 
their values at P, using superscripts (0) and (1) as above to denote orders in e and writing l a = + l^ a , 
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m a = m^ a + mM a . Using Eqs. (|T02|) . (|T03)) . (|T04|) . and (fT05)) . we have 

1 {0)T = ^, (106) 
/o 

^ (1)jf = 272 ( ? '- r o) + TT (*~*o) sin6 cos $, 

\/2r 

m« T = 0, 

72 1 

l (1)p = -o-(t-*o)- ^sin6cos$ + — (r-r ), 

,f2l n M (, 2 J 2 \ , \ 1 s JE(M — Tq) , s . n ^ M(rl + 2J 2 ), n2 
, JE(M - r ) (t _ Q $ _ 3A/_ (r _ ro)2 _ ro ^ Q + r_o ^ Q sin2 $ 



/ofoP 2r 2 / 2 p 2p p 

sin 2 6 cos 2 $ — "4" " (r — r )(t — t ) — "~ (r — r ) sin 9 cos $, 



£ 2 r .0^ 2^ 2MJ 2 . MJE 

— — sm 2 ecos 2 $- -j— r-r (t-t )- 77^- 



to (Dp = -^-(r 2 + J 2 cos<f>e- l *)sine-^^(t-t )e- 4 *, 
2r 0/ o 2r 0/ o 



m< 2 >" 



V^J^-^cos*. . . ^ V2JEM 



-/>!> 



(r-r ) sine 3 (t - io)(r - r ) 



2r 2 oP " " J " 2r 3 f p 
The terms in parentheses in Eq. (|53[) are then given to sub-leading order by 

2 



il T m"-m r V) 2 = (0 T m^-m^ T l^py 

+2 {l^ T m^P - m^ T l^P) (iWmW* - m^ T l^" + W - m« T Z«") ; (107) 



and p 5 is given by 

1 1 5 p( 3 ) 



p 5 p(2)f 2p(2)|' 



(108) 



implying 

3 m m„ fnYr,m„\ 2 . 15 m 



^° 2p(2)| V / 4 p( 2 )l V T J 

l(0)T m d)P _ m (0)T;d)p) ^(0)T m (2)p _ m (0)Tj(2)p + |(l)T m (l)p _ m (l)Tj(l)^ + Q^-ly (10Q) 



The expression for ipg as a mode sum is obtained from the value of this expression at t = to- The full expression, 
including terms involving t — to is given in Appendix ([C]). We denote by V , o _L an d 4'o~ SIj the leading and subleading 
parts of the singular field, respectively. Using Eqs. (|103[) . (|104l) and (|106|) . we obtain for to subleading order the 
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explicit form 



L 3m£ 2 r 2 sin 2 3m,/V 2,$ (r-r ) 2 3taE Je"** (r - r ) sin 

% - 72 + ^ f 2 P ' ^ ^ 

Jo P Jo r o P JO P 

_ SL 15mJ 2 Me- 2i * (r - r ) 5 15mJM.Be- 4 * sin0(r - r ; 

0o " ~ 
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15me- 21 * J 2 ((-i sin 4V 1 *) + ^ cos 2 $) (r _ ro) 3 gin 2 Q 

15mJSe- icI, (r 2 + 2 J 2 cos 2 $) (r - r ) 2 sin 3 8 15mr £ 2 (J 2 + r 2 + J 2 cos 2$) (r - r ) sin 4 9 
+ 1^ ^ + 2^ ^ 

9m (e -21 * J 2 M) (r - r ) 3 3me 4 * Jr E sin 3 6 15me- 4<I, A'/ (r - r ) 2 sin 6 
+ IFt W~ + To W~ + T 5 

9m J 2 (r-r )sin 2 e 

rofo P 5 

The axisymmctric part of each of these terms is proportional to an expression of the form 

S kl sin fc2 9 
( 5 2 + l-cos6) fc + 1 /2' 



(111) 



(112) 



where k\, ki and k are positive integers, and 5, given by Eq. (|C6[) . is proportional to r — rg. In Appendix [Cl following 
DMW, we use the generating function for Legendre polynomials, 

(e r + e J_ 2M)1/2 =E^ +1/2)|T1 ^)^ ^0, 

and its derivatives to express each term as a sum of Legendre polynomials and their derivatives. We then use a relation 
between the spin- weighted harmonics s Ye m and Legendre polynomials to write the series in terms of the harmonics 
2lfo- The leading order part of ipQ then has the form 



_ -m(r - 3M ) 3 / 2 / 1 \ 4tt« + 2)! 

( ^ o) - (0) - ri(r -2M)^ / h\ (*"2)!(2i+l) ° } ' (U4) 

where (0>q) is the axisymmetric part of 0q. 

oo 

Finally, each subleading term is proportional as a distribution to the sum ^2(1 + l/2)P;(cos0). That sum is a 

5-function with support at = [4g, and its projection along 2^0 therefore vanishes for all t. Eq. (|1 14[) thus gives 
(■05 ) to subleading order in L. 



C. Comparison with numerical determination of 0g 

We complete this section with a comparison of the analytic form of (ipo) with its numerical value obtained by 
matching to a power series in L. A comparison of the numerically and analytically determined contributions to 
the self- force from the axisymmetric part of tpo is given in Appendix [D] 

The retarded Weyl scalar tpffi is computed by integrating the radial Teukolsky equation for each value of £. The 
coefficients of 2^0 are matched to a power-series in L of the form 



4tt(£ + 2)! / B C 



^ = V (2* + l)(<-2)! [ A+ L + V + -)- < 115 > 

Shown below is a table of the fractional error in A and B when found numerically and compared to the analytic form 
given by Eq. (|1 14[) . Details of the numerical methods and checks of numerical accuracy are given in the companion 
paper. 
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ro/M 


-^analytic 


-^numerical 


|AA/A| 


|AB| 


8 


-0.002507548110466834 


-0.002507548108782573 


6 


717 


X 


10" 


10 


8 


496 


X 


10" 


10 


10 


-0.001214016915072354 


-0.001214016915092580 


1 


666 


X 


10" 


11 


1 


245 


X 


10" 


li 


15 


-0.0003356104323965837 


-0.0003356104323973954 


2 


419 


X 


10" 


12 


6 


003 


X 


10" 


13 


20 


-0.0001370231924969076 


-0.0001370231924969057 


1 


365 


X 


10" 


14 


3 


853 


X 


10" 


14 


25 


-0.00006882828667366571 


-0.00006882828667362643 


5 


706 


X 


10" 


13 


7 


462 


X 


10" 


15 


30 


-0.00003933557520091981 


-0.00003933557520088896 


7 


843 


X 


10" 


13 


1 


067 


X 


10" 


14 


35 


-0.00002455304484596332 


-0.00002455304484594233 


8 


549 


X 


10" 


13 


8 


341 


X 


10" 


15 


40 


-0.00001634080095822354 


-0.00001634080095821053 


7 


960 


X 


10" 


13 


5 


325 


X 


10" 


15 


45 


-0.00001141847437793787 


-0.00001141847437792919 


7 


605 


X 


10" 


13 


3 


673 


X 


10" 


15 


50 


-0.000008290448479296679 


-0.000008290448479290278 


7 


722 


X 


10" 


13 


2 


807 


X 


10" 


15 


55 


-0.000006208226467936966 


-0.000006208226467932644 


6 


961 


X 


10" 


13 


1 


912 


X 


10" 


15 


60 


-0.000004768831841073202 


-0.000004768831841069988 


6 


739 


X 


10" 


13 


1 


433 


X 


10" 


15 


70 


-0.000002990259098529844 


-0.000002990259098527988 


6 


209 


X 


10" 


13 


8 


488 


X 


10" 


16 


80 


-0.000001996831417921701 


-0.000001996831417920439 


6 


316 


X 


10" 


13 


5 


893 


X 


10" 


16 



TABLE II: The fractional error in the renormalization coefficient A and the error in B for (^>o) is given here for a particle in 
circular orbit in a Schwarzschild background at radius ro. A A and AB are the differences between the coefficients obtained 
numerically and by using the analytic expression ()C37[) . The analytic value of B is zero. 

V. BRIEF DISCUSSION 

The methods discussed in this paper for finding the self-force in a radiation gauge have been used to find the self-force 
on a particle in circular orbit in a Schwarzschild spacctimc, and work on orbits in a Kerr background is now underway. 
The advantage of a radiation gauge is the ease with which the retarded field can be computed. A disadvantage is 
the difficulty in analytically computing the singular field h s a p from ipQ. We have avoided this difficulty by using a 
numerical matching procedure to find the singuar field, and the companion paper shows that the numerical matching 
reproduces the renormalization coefficients for gauge-invariant quantities to machine accuracy, for the Schwarzschild 
example. 

We have shown that the perturbed metric in a radiation gauge generically has even parity to leading order in 
geodesic distance p to the particle trajectory. Using the renormalized field to compute the perturbed geodesic then 
relies on showing that the singular field gives no contribution to the self-force. The companion paper checks this for a 
particle in circular orbit in a Schwarzschild spacetime, in which the self-force is symmetric about a radial line through 
the particle. We numerically compute the axisymmctric part of the singular field and find that to sublcading order it 
coincides with the axisymmetric part of — mVp -1 ; it thus gives no contribution to the self-force at order p°. 

This regular behavior of the singular part of the self force may seem remarkable, given the line singularity in a 
radiation gauge that arises when one includes the perturbed mass in a radiation-gauge metric obtained from a Hertz 
potential. It is less surprising, however, if one considers a particle at rest in flat space. When principal null directions 
are chosen to have an origin not on the particle trajectory, ipg is nonzero, and the t > 2 part of the metric can be 
reconstructed in closed form in a radiation gauge using the CCK procedure. The self-force of course vanishes, with 
the contribution from the singular field coinciding with that from the retarded field; but the contribution from each 
is nonzero and gauge invariant under time- independent guage transformations. This implies that the singular part 
of the expression for the self-force in a radiation gauge has its Lorenz form -mV"/)" 1 - For a circular orbit, the 
result is implied by the fact that the gauge transformation of the self force can be written in a form, Eq. (|A24[) , that 
involves no derivatives of the gauge vector together with the fact that is O(p ). For generic orbits, the gauge 
transformation involves (u • V) 2 £ Q , and it remains to be seen whether the singular part of the self force can again be 
identified with — mV Q p -1 . 
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Appendix A: Gauge transformations 



1. Gauge transformation from Lorenz to radiation gauge 



We consider here the perturbed radiation-gauge metrics that describe asymptotically flat vacuum metric perturba- 
tions involving no linear change in mass or angular momentum. For harmonic time dependence, outgoing perturbations 
of this kind have, in a Lorenz (transverse-tracefree) gauge, the asymptotic behavior h a p = h a pe~ luJU 7 with 



and they therefore satisfy the IRG condition to 0(r~ 2 ): 

l^hap = 0(r- 2 ), / 



113 



0. 



0(r" 



(Al) 



(A2) 



We first find a corresponding asymptotically flat radiation-gauge metric perturbation by exhibiting a gauge trans- 
formation from the given Lorenz-gaugc metric to an asymptotically flat metric satisfying the exact IRG conditions. 

We then show that there is an asymptotically vanishing gauge transformation to an asymptotically flat ORG metric 
perturbation. Uniqueness of the gauge transformation implies that this coincides with the asymptotically flat ORG 
metric obtained from the ORG Hertz potential of Eq. (fH)]) . Here, for asymptotic flatness, we are requiring only that 
the tetrad components /i M „ have, for harmonic time dependence the form e~ tuu O(r~ 1 ); and for a time-independent 
perturbation involving no change in mass, angular momentum, or asymptotic dipole moment, = o(r~ l ). 
IRG metric for outgoing radiation. 

We obtain as follows the gauge transformation from a Lorenz gauge to an IRG metric perturbation, 



7 ret, RG i ret 



(A3) 



We begin with the transformation for a Schwarzschild background and then generalize it to Kerr. The transformation 
is described most simply in coordinates u, r, 6, 4>, with u the outgoing null coordinate. The harmonics of £ Q and h™p 
are then given by 



with the corresponding spin-weighted harmonics for the tetrad components of the metric perturbation, 



sYir, 



The radiation gauge condition l^h a p = has components 



Vi6 



■ V 2 fi 
V 3 £i 



-h IQt 
, n ll 



- -h rct 

— 'Ii2 

- -h rct 

— n 13 ■ 



These equations have the explicit form, 



(Or - -)|s 



1 



2 

w 



fd r 



M 



-hi2, 



l)/2] 1/2 |i 



with solution 



dr'huir') 
dr' 
dr' 



h^(r') + -f{r')h 11 {r')- 



£)6<o 



(A4) 
(A5) 

(A6a) 

(A6b) 
(A6c) 

(A7a) 
(A7b) 

(A7c) 

(A8a) 
(A8b) 
(A8c) 
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A result of Price et al. [36j shows, for a vacuum perturbation satisfying the condition l^h a p = 0, that the remaining 
gauge condition, h = 0, is also satisfied. That is, when hi^ = 0, we have h = — 2/134; and from Eq. (16) of that paper, 
the perturbed Einstein equation <5(Gn — 8ttTh) = implies 



i 34 = a 



/ e 



£ + * )+b°( e + g), 
q q; 



(A9) 



where a and b° are functions of u,9,<p. Now /134 = V3^4 + V4^3, and Eq. (|A8[) implies £3 = 0(r x ), whence 
^34 = 0(r~ 2 ). Since the right side of Eq. (|A9[) is 0(r~ 2 ) only if a = 6° = 0, we have /134 = 0. 

For each harmonic, it is not difficult to show that Eqs. (|A8[) give the unique solution to Eqs. (|A7I) for which 
the components hj^ vanish asymptotically for us 7^ and vanish faster than r^ 1 for u> = 0: Any other gauge 
transformation differs from the solution (|A8[) by a solution to the homogeneous equations, to Eqs. (|A7[) with = 0. 
Their general solution is 



£1 

fa 
fa 



-(iw + M/r)ki + k 2 , 
-[£(£+ l)/2] 1 ^k 1 + k 3 r. 



(AlOa) 
(AlOb) 
(AlOc) 



For uj nonzero, /122 vanishes asymptotically only if k± = and k 2 = 0; and /133 vanishes asymptotically only if k 3 = 0. 
Similarly, for ui = 0, /122 = o(r _1 ) only if fci = k 2 = 0; and /133 = o(r _1 ) only if fc 3 = 0. 

For a Kerr background, we similarly find the harmonics of the gauge vector in Kerr coordinates u, r, 0, cp of Eq. (|41l) . 
Harmonics of the gauge vector have the form 



£1 = fi(r, 0) e *(W-«™), £ 2 - f 2 (r, 0) e *("*-'"O, £ 3 - f 3 (r, 1 
The corresponding harmonics for the metric perturbation are 

The gauge transformation for a Kerr background is governed by the equations 



+ (A - 7 - 7)^ + (f - ^)£ 3 + (t - tt)U 
(6 - 2n)^ + (D + g)S; 3 



2 11 

-h ret 



= -h\ 



The components £ M are given successively by 



6 = - 



dr'hn 



dr g 



hta -(6- 27f)£i 



dr' 



ftl2 + (A - 7 - 7) fx + (iT - t)£ 3 + (7T — r)£ 4 



(All) 



(A12) 



(A13a) 

(A13b) 
(A13c) 



(A14a) 
(A14b) 
(Al4c) 



Asymptotic regularity follows from the asymptotic behavior (|Al[) of components along l a of outgoing waves in a 
Lorenz gauge. And the Price et al. result again implies h = 0. 



IRG perturbed metric for ingoing radiation 



The Hertz potential construction yields an asymptotically flat IRG form for each ingoing asymptotically flat metric 
perturbation. We find the gauge transformation from transformation from a Lorenz gauge to this asymptotically flat 
IRG. For simplicity, we restrict consideration to a Schwarzschild background. 
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In Eqs. (|A4I) and (|A5|) the outgoing null coordinate u is replaced by the ingoing null coordinate v = t + r* . In 



Eqs. (|A7[) <9 r is the replaced by 



. — 2imur* ft Jliuir 



-2iujr* 



d r e r , and the solution has the form 



dr'e 2iaJr /iu(r') 



6 = e 
fa = r 



-2iujr"' 



dr e 



h 1 2(r') + -f(r')h 11 (r')- + £ 



71/ 



^i»-^[^ + l)/2] 1/2 |i 



(A15a) 
(A15b) 
(A15c) 



Now, however, because the radiation is ingoing, h±x ~ e 1, 
relation 



D ikr* 



/r, when cj ^ 0. Asymptotic flatness follows from the 



+ 0(r~ 1 ) 



When = 0, asymptotic flatness follows from the asymptotic conditions (|A1[) . Again, because £3 = 0(r 1 ), we have 
V3^4 = 0(r~ 2 ), and the Price et al. relation then implies h = — 2/134 = 0. 



2. Gauge transformations of the self force 

A gauge transformation of the self- force was obtained by Barack and Ori [27|]- We give an alternate, covariant 
derivation, mention a second kind of gauge freedom, and obtain a simpler form of a gauge transformation of the 
self-force for a particle in circular orbit. 

A gauge transformation is an infinitesimal diffcomorphism that drags an unperturbed geodesic of the background 
metric to a neighboring curve that is a geodesic of the dragged-along metric. This can be stated precisely in terms of 
a congruence of timelike geodesies through a neighborhood of a point P. The dragged metric at P differs from the 
original metric by h a p = £^g a p, and the perturbed geodesic through P has 4- velocity altered by Su a = £^u a . The 
perturbed geodesic equation associated with a perturbation that is pure gauge has the form 

~8{u^Vpu a ) = £z(uPv p u a ) = 0. (A16) 

Writing 

£t(u p Vpu a ) = {£ i u fi )Vpu a + u' 3 V (3 £^ Q +u p [£ e ,Vp]u a (A17) 

and 

[£ € , V/jK = ^ 7 V^V 7 r - R a 7 p S u^ s , (A18) 

and with the perturbed acceleration defined by we have 

Sa a = -(u ■ V) 2 C + R a M5 u^v?i s . (A19) 

In this form, u a + 6u a is normalized to —1 with respect to the perturbed metric g a p + h a p, and the geodesic 
equation is affinely parameterized with respect to the perturbed metric. If the geodesic is parameterized so that its 
tangent is normalized to 1 with respect to the background metric, and we denote by 5^u a the change in u a with that 
normalization, then 

6 iU a = (Sf - u a up)8vP = (6% - u a up)£^. (A20) 
With <5 € a Q := S^V^u a - u^pS^u™, Eq. (|AT9]) implies 

S^ a a = - u a u f3 )(u ■ V) 2 ^ + R a p 7 6 U Pu' y Z 5 , (A21) 

and a a u a = 0. Note that the right side vanishes if £ a happens to drag a geodesic of the background spacetime to 
another geodesic of the background spacetime: This is the equation of geodesic deviation governing the connecting 



2G 



vector joining two neighboring geodesies of g a p. For general the right side of Eq. (|A21|) then measures the failure 
of £ Q to produce a geodesic of the background metric. 

The effect of a gauge transformation on the perturbed geodesic equation ((3|) is to replace 8u a by 8u a + S^u a and 
a a by a Q +(5 € a Q . 

There is a second kind of gauge freedom, an infinitesimal change in the background geodesic to which one compares 
a geodesic in the perturbed spacetime. The perturbed geodesic at an initial point of the trajectory can then be 
changed from u a to u a + Su a , with a a = 0. This allows one to regard the right side of Eq. (|A21[) as the change in the 
perturbed geodesic equation for a geodesic through an initial point P with the same initial tangent vector as that in 
the original gauge. 

For a particle in circular orbit, the gauge transformation of the self-force takes a simpler form for a gauge vector 
that is helically symmetric. Writing u a = u*/c Q , with k a the Killing vector t a + £l(f> a , and using the relations £k£ Q = 0, 
£k(& /3 V /3 £ Q ) = 0, and V (3 V 7 /c Q = R a ll3S k s , we obtain 

- (k • V) 2 ^ = i^V,3V Q (fc^ 7 ) - R a ^k^ s . (A22) 

The last term cancels the last term in Eq. (|A21[) to give 

S € a a = ^(^) 2 ^V^V Q (fcTfc 7 ), (A23) 

with corresponding gauge-transformed self-force 

.f Q = r + imCtWV/sV^fcT'fc,). (A24) 

For a particle in circular orbit in a Schwarzschild background, Eq. (|A23[) takes the form [4!| 

5 t a r = 3^7(1 - 3M/r)f . (A25) 



Appendix B: singularity and large i behavior 

Mode-sum renormalization involves relations, like the formal harmonic decomposition 



(1 



-1/2 



£=0 



-in 



1 + 1/2 



Y m {6), 



(Bl) 



between the short-distance (ultraviolet) singular behavior of a function and the large-^ behavior of its harmonics. The 
angle 9 is geodesic distance on the unit 2-sphere S to the origin 9 — 0, and in this example, a function that behaves like 
_1 has angular harmonics that behave like L -1 / 2 , where, as in the body of the paper, L = £ + 1/2. In this appendix, 
we review how one characterizes the large £ behavior of functions whose explicit angular harmonics are not known; 
and we give a precise meaning, in terms of distributions on the 2-sphere, to formal expressions like the divergent right 
side of Eq. (|Bip and to the angular harmonics of functions like / = 9~ n for which the integral J dClfY^ m diverges. 
We initially restrict the discussion to ordinary spherical harmonics and then generalize it to spin-weighted harmonics. 

The angular harmonics of the retarded fields are limits as r — > tq of angular harmonics of expressions that are 
nonsingular on spheres of radius r ^ r$. One therefore defines the angular harmonics of the singular field in the same 
way; and we end by showing that our formalism reproduces the harmonic decomposition defined in this way. 

The general relation between small-angle and largc-£ behavior is essentially identical to the relation between a short- 
distance singularity in a function /(x) and the behavior of its Fourier transform /(k) at large k. In one-dimension, for 
example, the n th derivative of a function /(x) is square integrable if and only if k n f(k) is square integrable, because 
/to - (ik) n f: 



I 



dx 



< oo 



dk 



k n f(k) 



< oo. 



(B2) 



Similarly, the n th derivative of a function f{Q,<f) = fimYim on S is square integrable if and only if l n fz m is 



square summable: With D a the covariant derivative operator of the metric ds = d9 + sin 



on S, the angular 
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Laplacian is D 2 := D a D a , and we have 

J dnD ai ■ ■ ■ D a jD a ^ ■ ■ ■ D a "f = J dnf(-D 2 ) n f < oc (B3) 

m+i)T\ft m \ 2 <™. (B4) 

One extends this relation to functions like (1 — cos^) -3 / 2 that are not square integrable by regarding them as 
distributions obtained by taking derivatives of functions like (1 — cosfl) 1 / 2 that are square integrable. If / is any 
distribution on S, fi m = J dQ,fY(, m exists, because Yi m is smooth. Thus, for example, writing 



{^D a D a + ^^)(1 - cos6»)t = (1 - cosfl)^ 1 (B5) 
or 2a 



gives / = (1 — cos6>) 3 / 2 as a distribution with 

f i0 = J (20(1 - cos6)- 1/2 (-2D a D a + l/2)Y m (already well defined) , 

or 

fio = J dn(l-cos6) 1 / 2 (2D a D a + 3/2)(-2D a D a + l/2)Y m . (B6) 

For any distribution /, the right side of Eq. (|B4[) is finite for some negative n, with increasingly singular distributions 
corresponding to increasingly negative values of n. To obtain a theorem that relates in this way the short-distance 
and large-^ behavior of distributions, we define the standard spaces of functions whose first n derivatives are square 
integrable, the Sobolev spaces H n , and then extend the definition to spaces H r of distributions, where r can be 
negative (and need not be an integer). 

Recall that the Hilbert space L2{S) is the completion of smooth (C°°) functions in the norm ||/|| defined by 



\f\\ 2 = J dO|/| 2 . The operator — D 2 + j is positive definite with eigenvalues L 2 . 



Definition. For positive integers n, the Sobolev space H n (S) is the completion of smooth functions in the norm 
defined by 



11/11* =Jdn 

Because the right-hand side is a sum of terms of the form 

C J dn.fD 2k f, k = Q,. 



(B7) 



a function has finite norm if and only if the function and its first n derivatives are square integrable. In particular, 
H (S)=L 2 (S). 

Because the operator — D 2 + j is positive definite, it has a well-defined square root, and we can write the definition 
of the norm in the more concise form 



1 



1/2 



\\f\\ n = \\V n f\\, with V := [-D 2 + . (B8) 

The action of the operator Dona distribution / is given by 

Vf = Y,Lf tm Y im . (B9) 

Equivalently, T>f is defined by its action on smooth functions g, 

Vf(g) = f(Dg), or f dQgVf = f dQ (Vg)f = £ Lg itn f im . (BIO) 
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We can now define H r (S) for all real r: 
Definition. The Sobolev space H r (S) is the completion of smooth functions in the norm || • || r defined by 



With inner product defined by 



the space H r is a Hilbert space. 
The relation 



ll/llr = |X> r /l ■ (BH) 

(g\f) r = [ dnV r gV r f, (B12) 



\m = J2 L2r \fe™\ 2 , (B13) 



Cm 



implies that a distribution / is in H r if and only if the sequence (fe m ) of its angular harmonics has finite norm 

\\(ft m )\\ 2 r :=J2 L2r \ft™\ 2 - (B14) 

Irn 

Formally, for each r, one can turn the set of sequences (fe m ) of complex numbers into a Hilbert space 4 with norm 
||(/em)||r- 

A function / is smooth if and only if it is in H r for all r, implying that the elements of a sequence (fe m ) are the 
angular harmonics of a smooth function if and only ]g m falls off faster than any power of I: Iim^oo i n \je,tr\ — 0- A 
second immediate consequence of the correspondence is the fact that T> maps H r to H r +\ and that the Laplacian 
maps H r to H r+ 2- (In fact, these maps are isomorphisms.) 

The function 9~ 1 is nearly square-integrable, with the integral J d£19~ 2 diverging only logarithmically and 
J d£l9~ 2+e finite for all e > 0. This suggests that 9~ 1 £ iJ_ e for all e > 0, and that is in fact the case: From 
Eq. (|Bip . the function / = (1 — cosfl) -1 / 2 satisfies 

ll(/£ m )i-e-^|^ m | 2 ^ 2e -47r^L- 1 - 2e <oo, alle>0. (B15) 

(Again, for e = 0, the sum diverges logarithmically.) Thus (1 — cos 6>)~ 1//2 G H_ t . Because [2(1 — cosfl)]" 1 / 2 differs 
from O^ 1 by O(0 2 ), the function 6~ 1 and any other function with the same singular behavior belong to f/_ e . From 
Eq. (|B5|) . successive applications of D 2 imply 9~ n G ff_„_ e . If a function has singular behavior 8~ n for integer n and 

if fim has singular behavior L s , for some s, it follows that s = n — —. 

Finally, a formal sum / = J2e m fomYi m (0, <t>), like the right side of Eq. dHTJ), has the meaning f(g) = ^Qtmfim, 
for smooth g. 

Finally we must relate this formalism to angular harmonics of functions f(r,9,(f>) that are singular at r = r and 
smooth for r r*o, when those harmonics are found as 

lim [ dnf{r,9,<t>)^ m . (B16) 

We suppose that / can be written in the form D r F, where F(r, 9, 4>) is continuous everywhere and smooth for 
r 7^ tq and where D is an operator for which D and D' have domains that include C°°(S). Then, for g smooth, 



J dflD n Fg = J dQFD^g =^> 
lim / dttD n Fg = lim / dttFD^g = [ dQ lim FD* 1 

•->r J r->r J J r-M-o 



9 



dnF(r )D^g(r Q ). 



4 The inner product is L 2r ge m fem> an d the Hilbert space is the completion in the norm ||-|| r of sequences for which limi^^ \fe m \i n = 
j 

(these are the sequences corresponding to smooth functions). 
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Because this last expression is, by definition, the action of the distribution /(ro) = D n F{ro) on g, we have the claimed 
equivalence 

hm f dSlfg = [ dtlf(r )g(r ). (B17) 



Appendix C: Mode-sum representation of (ipg) 

We present here details of the computation ipg as a mode sum, outlined in Sect. ([IV B[) . We first give the expression 
for i[>o to subleading order in Schwarzschild coordinates, including terms involving t — to. As in Sect. (|IV B[) we denote 
by 4>o by ipQ L and ipQ Sh , respectively. 

Eqs. ()103|) . (|104|) and (|106|) yield for to subleading order the explicit form 



3me 



-2i* 



/oSV 
15m 



(£r 2 sin9e I$ + J(r - ro) - J/o(t - Me" 4 *) 



4p 7 
2ME 2 



( r ~ r o) 7i2 sin0 O - r o) + 7T~i (* - *o)(r - r ) 4 



J-4 4 
JO 

(r ~ r ) 3 sin 2 9 + 



f4 2 

Jo'o 

2J 2 (2J 2 cos 2 $ + r 2 )e- 2 



f2 3 

Jo'o 



fo Jo'o 
4JEe" 2i *(e i *Mr 2 + J 2 (A/ - r ) cos$) 



f a r 4 
jo'o 

-(r-r ) 3 sin 2 9 + 



4e 



-2i$ 74 



J 4 M 



J-2 6 

Jo'o 



(r-rof(t-t )' 2 



f3 4 
JO'O 



s . ^ 4J£e- l *(r 2 + 2J 2 cos 2 $) / x2 . 
(r - r ) 3 (t - t ) sin 9 + ^ -(r - r ) 2 sin 3 9 



4 j2 e -2i$[ ro | r() _ M + y or() + e 2 1 *(r„-A/)|^2 + j2 j2 CQf , 2$ ] 
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2e- j;,I, J 3 75;((5 
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ro/o 
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f 2 r 3 
Jo'o 



4e 
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(CI) 
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We repeat the form at t = to given in the text, in order to label each term with a subscript for later reference. 
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15m JMEe' 1 * sin9(r - r ) 4 



P 5 



(C2) 
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(C3) 



The axisymmetric part of each of these terms is to be written as a sum over 2^ff) (®> 0) at r = r . Each term in the 
singular field involves the leading part of p 2 , namely 



where 



.4 



A(r - r ) 2 + B(l - cos 6) = \ t=ta + 0(9 4 ), 



B = 2r, 



r - 2M 



r - 3M 



r - 2M 1 

We follow the notation of DMW, writing 

p 2 := B{6 2 + 1 - cos8), 5 2 := A(r - r ) 2 /B. 
Then the leading term in ipQ is given by 



M sin 2 $ 
r - 2M ' 
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(C5) 



(C6) 



(C7) 



The axisymmetric part of the above expression is achieved by angle averaging over Substituting the values of A, 
B, E and J in the above expression, we obtain 



(^0" L ),o(©) 
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S sin 9 



(C8) 



where (/($)> = (2*)- 1 /(*)#&• 



We start with a form of generating function of the Legendre polynomials given by 

1 



£ e -(m/2)m W) TjtQt 



(C9) 



(e T + e- T - 2u)V: 

and set it = cos 9, |T| = a/25. In the limit T ^0 ± , the sum does not converge, but it is well-defined as a distribution 



tTo (e T + e- T - 2«)Va 



(CIO) 
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where the symbol = means equality of both sides as distributions on the sphere. In particular, the regularization 
proceeds by imposing a cutoff £ max on the singular field and on the retarded field; the projection V of the distribution 
(|C9j) onto the subspace I < £ max is the smooth function 



li m r 1 

T^o (e T + e- T - 2w) 1 /2 



E p *( u )> 



(Cll) 



z=o 



and the subsequent functions obtained by taking derivatives can be regarded as projections of the corresponding 
derivatives of the distribution (|C10|) . In particular, taking successive derivatives with respect to T gives the relation 
(Eq. (Dll) of DMW) 



lim 



t^o (e T + e~ T - 2«) fc +V 



1/2 E 



(21+1) 



2(2fc- l)T 2k - 



(C12) 



sin" 



Consider now the expression (|C7[) for the leading term in i/jq. The first term is proportional to ^2 +1 _ cos q)b/2 ; 
where S is proportional to T. We express this term as a sum of Legendre polynomials by differentiating (|C9[) twice 
with respect to u to obtain 



1 



(e T + e- T - 2uf/ 2 



1 OO 

l£e-<*-wjy( u ). 



(C13) 



i=o 



The second term is proportional to 



(5 2 + l-cose) 5 /2 



The fact that it is 0(e 3 ) suggests that it can be written as a 



linear combination of derivatives of -r^r-, — ^ — /, , 



oid u 



1 



1 



(e T + e- T - 2w) 5 / 2 """(e T + e- T -2w) 1 /2 1 ^ T ( e T + e -T _ 2tt )i/2 ' 
and solving for a and (3 gives a = /3 = jxj- Then 

T 2 



(e T + e- T - 2m) 5 /2 2 + 



4- 7/ ^-^ 



-(i+l/2)T 



(C14) 



(C15) 



lim 



2 3/2 



E 



pi 



s^o (5 2 + 1 -cos6) 5 / 2 3 
The axisymmetric part of the third term (its angle average over $) vanishes 



(C16) 



We next use the same techniques to express the subleading terms in x/Jq as power series in Pi with coefficients 
polynomial in £. The axisymmetric parts of the second, fourth, seventh and eight terms vanish. In each of the 
remaining terms, Eq. (IC12I) is used to expand an expression involving a power of S 2 + 1 — cos 0, and in each case we 

OO 

find that the term is proportional as a distribution to the sum 

EG + l/2)Pz(cos6>); that sum is a (5-function with 

support at 6 = |4|: 



EG + l/2)Pi{cos0) = 5(1 - cos( 



(C17) 



Because 2^,0 vanishes at = 0, the expansion of the singular field in terms of 2^0 has n o subleading contribution. 

We verify the claimed form of each of the subleading terms in Eq. (|C3|) as follows: 
The first term is proportional to <5 5 /(<5 2 + 1 — cos0) 7 / 2 , and Eq. (|C12|) gives 



S 5 2 / 1 \ 

lim — — = -Vh?+- /Vcos0). 

<5^o+ (J 2 + l-cos0) 7 /2 5^V 2/ 



(C18) 
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The third term is proportional to 

<5 3 sin 2 9 



2<5 3 



2S 5 



(<5 2 + l-COs9) 7 / 2 (52 + 1 _ COS 6) 7 /2 (52 + 1 _ cos 9)7/2 + ( £ )> 



and Eq. (|C12|l gives 



lim 



S 3 sin 2 9 



«5_,o+ (<5 2 + l-cos6) 7 / 2 15 
The fifth term is proportional to 

8 sin 4 9 4(5 



4<5 5 



4(5 3 sin 2 9 



(5 2 + i -cos e) 7 /2 (J2 + i -cos 6) 3 / 2 (<5 2 + 1 -cos6) 7 / 2 (<5 2 + l-cos9) 7 / 2 ' 
and Eq. (|C12|) gives 

<5sin 4 9 



lim 



G4 



^0+ (S 2 + 1 - cos6) 7 / 2 15 



£ 



1 



Pf(cos9). 



The sixth term is proportional to 



lim 



S-S-0+ (S 2 + 1 - cos9) 5 / 2 3 
The ninth term is proportional to 

<5 sin 2 6 2<5[l-cos9 + 0(9 4 )] 25 



,£(W 9 Wos9) 



2c5 3 



(52 + 1 -cos e) 5 / 2 ((5 2 + l -cos e) 5 / 2 (5 2 + 1 - CO s9) 3 /2 (J2 + 1 -COS 9)5/2 
again using Eq. (|C12| we have 



5sin 2 9 . 8^/ 1\ . ^. 

I™ 77^ -7777 = - > K + - iVcos9). 

5^o+ (<5 2 + 1 _ cos 9 5/2 3 ^ V 2/ 



(C19) 



(C20) 



(C21) 



(C22) 



(C23) 



(C24) 



(C25) 



The axisymmetric parts of the second, fourth, seventh and the eighth terms vanish. 



We have obtained the leading terms as series of Legendre polynomials, and we now convert them to series involving 
2Y1Q. We begin with Pj . We use the relation between s Y^ m and D e _ sm (9, $, 0), the representation matrix for the 

rotation group [50j to write 2^0 = Yn- We next convert the series in terms of Legendre polynomials to series involving 

1/2 

Yim-, we have 



2 Y m . Then, using PV 



Att (l+m)\ 
(21+1) (l-m)\ 



pf ) (cos9) =Y,Cw 2^'o(9,0), 



where 

_ Un(£ -l)l(£+ l)(£ + 2y 1/2 
~ [ (2£+l) 

Regarding P^(cos9) and P f '(cos9) as elements of L2(S 2 ), we find 5 



P £ ( C 0s9) = Y, A ^ P nHcOSe), 

n=2 
00 

P;(cos9) = / T]P f „P I l 2) (cos9), 



(C26) 

(C27) 

(C28) 
(C29) 



Tl=2 



5 That is, the coefficients Ai n and Bg n are obtained from the inner products of with P^(cosO) and P^'(cosQ), and the relations are 
implied by L2 completeness of iYt m . 
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where 



A( n 
Bfn 



(2n + l)(n-2)! 

2(n + 2)! 
(2n + l)(n-2)! 

2(n + 2)! 



2^-1) 
' (2£+l) t? 

2^ + l)d,_i, n , 



(C30) 
(C31) 



with 



1, n-f a positive even integer, 
0, otherwise. 



The form of these coefficients was found by numerical experiment for £, n < 10 and then checked as rational numbers 
for larger values. 

The second term in the expression for the leading part of the singular field is proportional to the right side of 
Eq. (|C16|) ; and we now show for each £ that the bracketed expression vanishes as an element of L 2 - We have 



'=0 L 



C«' 2^'0- 



(C32) 



For £' even, the sum over £ in Eq. (|C16|) is proportional to 



2{£' + 2)! 



(2£' + !){£' -2)! 



E 



Apt, 



B/fi 



The second and third sums are given by 

l'-2 



e=o 



e'/2~i 



£=0, even 



E 2^+1) 



1=1. odd 



E (2^+l) 2 = £ (4p+l) 2 = V(4£ ,2 -6£'-l), 



(C33) 



(C34) 



^—0, even 



p=0 



«=1, odd 



£ 2£(£+l) = [2(2p+l) 2 + 2(2p+l)] = + 2)(2£' - 1), 

p=0 



(C35) 



and the identity 



E 

1=0 



An> + \£+-\ Bw 



0, for £' even, 



(C36) 



follows. A similar manipulation yields the same identity for £' odd. We conclude that the projection of the distribution 



E 



=0 L 



Pi + I I + g I ^ 



along 2 Yeo vanishes. 

Then, of the three terms in Eq. (|C7[) for (ipf,), only the term proportional to sin 2 0//5 5 / 2 is nonzero when written 
as a sum over 2^0; and the axisymmetric part of V'o nas j to subleading order, the form 



-TO 



(r - 3A/) 3 / 2 / 1 



r 2 (r -2M)5/2 \ % 5/ 



E 

1=2 



' 4tt(£ + 2)! 
(^-2)!(2^+l) 



2 ^o(e,o). 



(C37) 
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Appendix D: Comparison of analytic and numerical computation of self-force 

To show how accurately one can recover the leading and the subleading terms in L in the mode sum expression for 
a s r by numerically matching a power series in L to the values of a£ et , we will present an example where we know A 
and B analytically: the contribution to the self- force from the part of h\\ that is axisymmctric about a radial line 
through a particle in circular orbit in a Schwarzschild background. The contribution to the self-acceleration from hn 
is 



i[hu}) = 



(1 - 2M/r) 2 
8(1 - ZM/r) 



2A 



1 



2M 
r 



d r — 6 



M 



(hi 



(Dl) 



To compute the leading and subleading terms analytically requires us to find the leading and subleading terms of 
the radial and time derivatives of tpQ. From Eq. (|C1|) , we have 



($Vo (* = *»)> = 



/ 6^J 2 (r-ro)e- 2 ' $ 
IhpJE sin 9 cos <f> 

W 1 



6/x£;Jsin 3 ee 



hp 



r- 



E 2 r^ sin 2 9 + — (r-r ) 2 e 



2 -2i$ 



2JE sin9(r - r )e 



— a 



2x7 



(D2) 



with only two terms surviving the angle- average over <£>. The first term is proportional to ^ S 2 +1 S _ u ^/2 where 5 cx (?* — r ) 
and u = cos 9. Taking d u &r of Eq. (|C12[) . we have 



lim = — 

*_*)+ (5 2 + 1 - uf/ 2 3 



^(£+l/2)P;( u ) = (D3) 



The second term is proportional to ^ S 2 + S 1 s ^. Q ^ e y/2 ■ A straightforward calculation gives us 

Bm Jsin2Q _ =-Y(l + l/2)R (2) (cos 9) (D4) 

8^o+ (<5 2 + l-cos9) 7 /2 15 ^ v 7 ; f v ; v ; 

Therefore, we get 



X ^V — (T^i — 2 ^ o(0) 

A similar calculation for the radial derivative of the leading singular field gives us 



^(r - 3A/) 2 / 1 \ A*M(r - 3M) 



-2i* 



v / 4^(21 +!)(£ + 2)! 

Here the superscript S-L refers to the singular leading. The dot and prime represent time and radial derivatives, 
respectively. 

To find the t and r derivatives of the subleading terms, we need to assume the following result, directly verified 
only for small n. 

6 2n 

lim — r = 0, n > 1. (D7) 

5 ^o+ (S 2 + 1 -cos9)"+ 3 / 2 ' v ' 
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Using Eqs. (|D2|) and (|D3|) . and (|D7|) , we have 

f5fiM(r - 3M) 3 / 2 /cos 2 $\ 3^M(r - 3Af) 3 / 2 / 



V 4(r a -2M) 5 /2 \ ? 

' 4tt(^ + 2)! 
(2^+1)0? -2)! 



7/2 / r4(r -2M) 5 / 2 \x 5 / 2 

2 ^,o(e) 



(D8) 



SL / 15 M (r - 3M) 5 / 2 //sin 2 3>\ r - M /cos 2 $\\ 3^/(r - 3M) 3 / 2 / 1 \\ 

Wo ) t0 ,r^r \2V/ 2 r%{r -2My/ 2 \\ X 7/2 / r -MI\ x 7/2 / J r*(r - 2AI)V 2 \x 5/2 // 
v - / 4tt(£ + 2)! " 

Here the subscript S-SL refers to singular subleading. 
Prom Eqs. ([iO]). ([32]) and (|Dl) . we find 



with 



(o[/in]) = AL + 5 + 0(ir 2 ), (D10) 



^ = -5ciM + 4(4c 2 M + c 5 (r - 3M)) 
8ro /2 (r -2M)V2 



where 



- 2 i<l> 



(r - 3M) 1 / 2 (-20c 3 M + 2c 4 (7M - 2r ) + 5c 6 (r - 3M)) 

, - 2A/)V2 

1 \ _ / 1 \ _ / sin 2 g \ , t-q-M / cos 2 $ 



S - 2r 3 (r -2M)V2 



In using Eq. (|4T)]) to calculate A and £?, we ignore the term involving oc riifl^S, because it smaller than the leading 
term by three powers of I. 

We numerically calculate (a r,Tet [hn]) by matching it to a scries in L of the form 

^T,o= AL + B + ^ + § + ---. (D13) 

Shown below is a table of the fractional error in A and B when found numerically for the self-force's contribution 
from the axisymmetric part of hi±. 
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Af 


-^analytic 


-^numerical 


|AA/A| 


^analytic 


-^numerical 


\AB/B\ 


8 


0.01025269710132717 


0.010252697099652308 


1.634 


X 


10" 


io 


0.005653763844987715 


0.005653764469505427 


1.105 x 10" 


-7 


10 


0.006104203254814939 
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X 


10" 


11 


0.003815279509214152 
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TABLE III: The table compares values of regularization parameters calculated analytically to values obtained numerically by 
matching the retarded field to a series in (£+ 1/2); the quantity that is regularized is (a r ' ret [/ui]), as described in the text. The 
first column lists orbital radius in units of Schwarzschild mass; the second and the fifth columns list the analytically computed 
leading and the subleading regularization parameters A and B; the third and the sixth columns list the numerically obtained 
values of A and B, and the fourth and the seventh columns list fractional differences between the analytic and numerical values. 
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